Hamiltonian thermodynamics of charged three-dimensional dilatonic black holes 
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. The action for a class of three-dimensional dilaton-gravity theories, with an electromagnetic 

Maxwell field and a cosmological constant, can be recast in a Brans-Dicke-Maxwell type action, 
with its free to parameter. For a negative cosmological constant, these theories have static, electri- 
cally charged, and spherically symmetric black hole solutions. Those theories with well formulated 
asymptotics are studied through a Hamiltonian formalism, and their thermodynamical properties 
are found out. The theories studied are general relativity (to — > ±oo), a dimensionally reduced 
cylindrical four-dimensional general relativity theory (w = 0), and a theory representing a class of 
theories (to = —3), all with a Maxwell term. The Hamiltonian formalism is setup in three dimensions 
through foliations on the right region of the Carter-Penrose diagram, with the bifurcation 1-sphere 
as the left boundary, and anti-de Sitter infinity as the right boundary. The metric functions on the 
foliated hypersurfaces and the radial component of the vector potential one-form are the canonical 
»' ! ' coordinates. The Hamiltonian action is written, the Hamiltonian being a sum of constraints. One 

finds a new action which yields an unconstrained theory with two pairs of canonical coordinates 
{M,Pm; Q,Pq}, where M is the mass parameter, which for lo < — | and for ui = ±oo needs a 
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careful renormalization, Pm is the conjugate momenta of M, Q is the charge parameter, and Pq is 
its conjugate momentum. The resulting Hamiltonian is a sum of boundary terms only. A quanti- 
zation of the theory is performed. The Schrodinger evolution operator is constructed, the trace is 
taken, and the partition function of the grand canonical ensemble is obtained, where the chemical 
potential is the scalar electric potential <f>. Like the uncharged cases studied previously, the charged 
black hole entropies differ, in general, from the usual quarter of the horizon area due to the dilaton. 
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I. INTRODUCTION 

In a previous paper flj we have motivated and studied the Hamiltonian thermodynamics of three-dimensional 
dilatonic black holes. To construct a classical Hamiltonian formalism is important for several reasons. First, it is always 
elegant to write the equations of the theory as a pair of symmetric first order equations. Second, if one can put the 
classical theory into a Hamiltonian form, then by applying certain rules one can get a quantized version of the theory in 
a first approximation. Third, through Euclideanization of the Schrodinger time evolution operator exp(— iHt), where 
H is the Hamiltonian, one gets the partition function, which in turn leads to a thermodynamical decription of the 
system. Thus, the Hamiltonian formalism and the thermodynamical description are linked subjects. The Hamiltonian 
thermodynamics of several different black hole systems have been studied first by Louko and Whiting in Q and further 
analyzed in @, H, H, H, 0| and [H, (see also The Louko- Whiting method relies on the Hamiltonian approach of 

[Io| . which in turn is an important ramification of the Arnowitt, Deser, and Misner approach, the ADM approach, and 
its major developments [Tl|,[l3|. Other methods, like direct calculation from quantum fields in curved spacetime (see, 
e-g, or path integral methods (see, e.g, [H,[lB, E3, 0, OH]), have been used to study the thermodynamics of 
black holes. Now, the study of thermodynamics of black holes in any specific dimension is important, for instance, to 
understand universal properties independent of the dimension itself. In particular, one can single out three dimensions 
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as an adequate interesting dimension, because three dimensional general relativity has a special black hole with simple 
properties [20], and dilatonic extensions of this theory have black holes with properties similar to four-dimensional 
general relativistic black holes [U, ■ 

The main purpose of this paper is to generalize our previous work on Hamiltonian thermodynamics of three- 
dimensional dilatonic black holes [H, to black holes, still in three dimensions, in which the graviton-dilaton Brans- 
Dicke theory, with its uj parameter, is coupled to a gauge field, namely, the electromagnetic field. Three-dimensional 
Einstein-Maxwell theory with its black hole solution has been studied in [13] (see also [20]), and electric dilatonic 
extensions also have charged black holes (23. [25| (see also (2lL l22l|). Further analyses on the properties of these black 
holes have been worked out in [26l. I27L l28| . In principle, ui can take any value, i.e., oo > u> > — oo. But, since we want 
to perform a canonical Hamiltonian analysis, using an ADM formalism supplied with proper boundary conditions, it 
is necessary to pick up only those ujs which yield black hole solutions that fulfill boundary conditions adequate to a 
Hamiltonain method. As in the uncharged case the cases of most interest are then black holes for which uj — ► ±oo, 
uj = 0, and ui — —3. Now, unlike the uncharged case, it was realized in (23l . [25| (see also [2(j) that the inclusion 
of electric charge creates a complication to the definition of the mass M. Indeed, for uj < — §, which includes the 
general relativity case (since it can be considered as the uj — > — oo limit), one finds a divergent asymptotic behavior 
of the mass. This problem can be cured by a proper mass renormalization (23l . [25| . But the whole issue has never 
been confronted when one makes a full use of the Hamiltonian formalism, as we do here. We solve this difficulty and 
show that the mass redefinition conforms wholly with the formalism for the ui — > oo case, is not needed in the uj = 
case and it is compatible with the ui — — 3 case. In each case, the black hole system and its thermodynamics can then 
be studied. We have chosen to study the systems in a grand canonical ensemble. In such an ensemble, the boundary 
radius I, the temperature T (or its inverse (3), and the electric potential (i.e., the chemical potential) 4> at the boundary 
radius are fixed. The electric potential is conjugate to the electric charge and, in this ensemble, the charge itself is not 
fixed. Our main results are, a complete derivation from a Hamiltonian formalism of the mass-energy (E), the entropy 
S, and the charge (Q) , as well as the establishment that the grand canonical ensemble for charged three dimensional 
black holes is locally stable for the black holes studied. Thus, the choice of such an ensemble is appropriate, since the 
corresponding boundary onditions lead to a well posed problem and a stable ensemble. 

The structure of the paper is as follows. In Sec. |ll]we present the classical solutions of the three-dimensional dilatonic 
charged black holes, whose quantization through Hamiltonian methods we will perform. There is a free parameter 
uj for which we choose three different values, uj — oo, 0, —3, corresponding to the BTZ black hole, the dimensionally 
reduce four-dimensional cylindrical black hole, and a three-dimensional dilatonic black hole, respectively. We also 
introduce the spacetime foliation through which we will define the canonical coordinates and which will allow us 
to write the action as a sum of constraints multiplied by their respective Lagrange multipliers. Then, follow three 
sections, Sees. IIIII IIVI and|Vl where we develop the thermodynamic Hamiltonian formalism for the uj = oo, 0, —3 
charged black holes, respectively. We divide each section into A-F subsections to better analyze the Hamiltonian 
thermodynamics of the systems. 



II. CHARGED BLACK HOLE SOLUTIONS ALLOWING HAMILTONIAN DESCRIPTION 

A. The 3D black hole solutions 

Three-dimensional theories which contain black holes have been studied in several works. Indeed in (20l. l2ll. I22L I27I 
HH, a Brans-Dicke action, with gravitational and dilaton fields plus a cosmological constant, was used to study such 
solutions. Now, we can couple the graviton and dilaton fields to an electromagnetic field, in order to find charged 
black hole solutions in three dimensions. For three-dimensional general relativity this has been done in (see also 
[20lp. whereas for Brans-Dicke theory this has been done in [25] (see also [28]). The Brans-Dicke action with a Maxwell 
term is given by 

S =^J d 3 x^e~ 2 ^ (R - Md<f) 2 + 4A 2 + F^F^) + B , (1) 

where g is the determinant of the three-dimensional metric g M „, R is the curvature scalar, <j> is a scalar dilaton field, 
ui is the Brans-Dicke parameter, A is the cosmological constant, F^ v — d^A v — d v A^ is the Maxwell tensor, where 
A = A^dx 11 is the vector potential one-form, and B is a generic surface term. A general solution is a solution for the 
metric, the dilaton and the gauge field. We search for static solutions. The solution for the metric is then [HI 
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+ ^\nR)R 2 -bR dT 2 + -, r + R 2 dtp 2 w = -2 , (3) 

4 / J M + ^lnR) R 2 -bR 

*> - (uw tocuo + *"*•) - 4A3g ,„g ) + x2R4 + *V. — 1 , w 

ds 2 = -dT 2 + dR 2 + dtp 2 , w = -l, (5) 

where T, i? are Schwarzschild coordinates, a is a constant related to the cosmological constant (see below), and 6 is a 
constant of integration related to the mass (see below), x ls constant of integration related to the electric charge, and 

k is defined as k = f ^+2) • ^ ne g enera l solution for (j) is given by 

= —L^lnM), ^-1, (6) 

= constant, u = —1 . (7) 
The general solution for the vector potential A, written as A = A^{r)dx^ — At(t) dT, is given by 

A = ^x(u + l)(aR)-^dT, (8) 

A = 0, lu = -1. (9) 



For the constant a one has 



2|(W + 1)A| ^-2,4-1. UO) 



v /|( W + 2)(2c + 3)| ' ^ ' 2' 

a = 1, c^ = -2,-|, (11) 

a = 0, w = -l. (12) 

For a; = —2 we have A 2 = — x 2 /16, which means that, unlike for the uncharged case, the cosmological constant is not 
null. The constant b is related to the ADM mass of the solutions by 

M = ^b, -|, -1, (13) 

U) + 1 2 

M = -b w = -2, (14) 

M = -4A 2 lno, w = ~\i ( 15 ) 
M = 0, lj = -1. (16) 

The constant x is s related to the ADM charge of the solutions by 

X 2 = Q 2 , c^-2,-|-l, (17) 
X 2 = 0, u = -l. (18) 

For w = —2,-1, there is no good definition for Q in terms of x- 

Since we want to perform a canonical Hamiltonian analysis, using an ADM formalism supplied with proper boundary 
conditions, it is necessary to pick up only those solutions that fulfill the conditions we want to impose. First, we 
are interested only in solutions with horizons, so we take b to be positive. Second, apart from a measure zero of 
solutions, all solutions have a non-zero |A|. This does not mean straight away that the solutions are asymptotically 
anti-de Sitter. Some have one typ e or another of singularities at infinity, which do not allow an imposition of proper 
boundary conditions. So, from [25J with the corresponding Carter-Penrose diagrams we discard the following solutions: 
u) = — 1 which is simply the Minkowski solution, — 1 > u > — | since it has a curvature singularity at R = +oo, which 
is inside the horizon, and has a null topological singularity at R = 0, cj = — 1 since all the Carter-Penrose boundary 
is singular, the same happening in the interval — | > u> > —2. Thus, the cases of interest to be studied are black holes 
for which oo > lu > — 1 and — 2 > uj > — oo. For b positive these solutions have ADM mass M positive. As in []]] we 
choose three typical amenable cases where an analytical study can be done. The first case is the charged BTZ black 
hole [13] (see also [20]), which can be found by taking appropriately the limit ui — > ±oo. This black hole is thus a 
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solution of the Einstein-Maxwell action in three dimensions. The other cases are w = and u> = —3. The theory with 
lu = is equivalent to cylindrical four-dimensional Einstein-Maxwell general relativity, and the corresponding charge 
black hole was found in see also (2j|. The theory with lj = — 3 is just a case of 3D charged Brans-Dicke theory, 
with a black hole solution that can be analyzed in this context [25]. 

Now, as mentioned in the introduction, the definition of the ADM mass M is not straightforward when there is 
electric charge, in the sense that one can face a divergent asymptotic behavior of the mass M, see [H, H3| as well as 
[2o| . The problem is more acute when oj < — |, which includes the gen eral relativity case (since general relativity can 
be considered the u> — > — oo), see [25| for a full discussion (see also |2Cl l23|). For such oj, the mass M can be written 
as M = Mq = q + DivM(Xj-R)> where Mq = q is the mass for the uncharged case and ~D'ivm(Xi R) ls the divergent part 
when there is charge. To treat this divergence, as explained in [25], we consider a large radius i?* frontier. Then we 
write the mass M as 

M = M(R*) + [Div M (x, R) ~ Div M (x, #*)] , (19) 
where the function M(R*) is equal to 

M(R t ) = M Q=0 + Div M (x,#*). (20) 

The term in square brackets in Eq. (fl9|) tends to zero as the radius R tends to the frontier radius, R — > -R*. M(R^) 
is interpreted as the energy inside the radius i?*. From Eq. (|20|) the term — DivM (Xj R*) ls interpreted as the energy 
outside the frontier in i?*, apart from an infinite constant, which is hidden in M{R ie ). The total mass definition, M, 
is then well defined. In practical terms, one needs not pay attention to the divergent factors DivM because all one 
needs to do is consider null the asymptotic limits, limDivM(x, R) = 0, ignoring the frontier at i?*. For oj > —3/2 
the mass term behavior of the three-dimensional black holes is similar to the Reissner-Nordstrom case, and it offers 
no problems. Indeed, for the Reissner-Nordstrom black hole the metric function gu and the mass term, following the 

prescription above, are written as gu = 1 — 2M jfi - + Q 2 — -j^) , and M — M(R* ) + respectively. The energy 

o 2 

inside the frontier i?* is then given by M(R*), and the electrostatic energy outside is given by As this last term 
tends to zero as i?* — > oo there are no real divergences in this case, thus implying that M(R*) carries no infinity 
within, and so the prescription is trivial, it gves nothing new. This is similar to what happens for three-dimensional 
black holes with oj > —3/2. On the other hand, for u) < —3/2 one should follow the prescription carefully, i.e., one 
has to hide an infinite constant in the definition of the energy within the frontier at i?* . 

B. ADM form of the metric and dilaton 

The ansatz for the metric and dilaton fields with which we start our canonical analysis is given by 

ds 2 = -N(t, rfdt 2 + A(i, r) 2 {dr + N r (t, r)dtf + R(t, r) 2 dip 2 , (21) 



^ = (aR(t,r))^ . (22) 

where {T, i?, ip} are the Schwarzschild type time, radius, and angular coordinates. This is the ADM ansatz for 
spherically symmetric solutions applied to three dimensions. In this we follow the basic formalism developed by 
Kuchaf The canonical coordinates R and A are functions of t and r, i.e., R = R(t,r) and A = A(t,r). Now, 

r = is generically on the horizon as analyzed in (Kjjbut for our purposes r — represents the horizon bifurcation 
point of the Carter-Penrose diagram [3] (see also [3|-@l an d [H)- In three spacetime dimensions the point actually 
represents a circle. Asymptotically the coordinate r tends to oo, and t is another time coordinate. The remaining 
functions are the lapse function N — N(t,r), and the shift function N r = N r (t,r), and they will play the role of 
Lagrange multipliers of the Hamiltonian of the theory. The canonical coordinates R = R(t,r), A = A(t, r) and the 
lapse function N = N(t, r) are taken to be positive. The angular coordinates are left untouched, due to spherical 
symmetry. The dilaton is a simple function of the radial canonical coordinate, and it can be traded directly by it 
through equation l(22|) . as will be done below. The ansatz f2"Tj) - ([22"|) is written in order to perform the foliation of 
spacetime into spacelike hypersurfaces, and thus separates the spatial part of the spacetime from the temporal part. 
Indeed, the canonical analysis requires the explicit separation of the time coordinate from the other space coordinates, 
and so in all expressions time is treated separately from the other coordinates. It breaks explicit, but not implicit, 
covariance of the theory, but it is necessary in order to perform the Hamiltonian analysis. The metric coefficients 
of the induced metric on the hypersurfaces become the canonical variables, and the momenta are determined in the 
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usual way, by replacing the time derivatives of the canonical variables, the velocities. Then, using the Hamiltonian, 
one builds a time evolution operator to construct an appropriate thermodynamic ensemble for the geometries of a 
quantum theory of gravity. Assuming that a quantum theory only makes sense if its classical form can be quantized 
by Hamiltonian methods, one should pick up only solutions which can be put consistently in a Hamiltonian form. 
Thus, in the following we perform a Hamiltonian analysis to extract the entropy and other thermodynamic properties 
in the three-dimensional Brans-Dicke-Maxwell black holes mentioned above, those for which u — oo, 0, —3. 



C. ADM form for the vector potential A 

The canonical description of the vector potential one-form is given by 

A = Tdr + ^dt, (23) 

where V and $ are functions of t and r, i. e., T(t, r) and <I>(t, r). The function T(t, r) is the canonical coordinate 
associated with the electric field, and the function $(t, r) is the Lagrange multiplier associated with the electromagnetic 
constraint, which is Gauss' Law. In the same way as with the gravitational degrees of freedom above, when the Maxwell 
term in the action (U) is written, the conjugate momentum to the coordinate Y(t, r) can be derived the usual way 
from the Lagrangian. 



III. HAMILTONIAN THERMODYNAMICS OF THE CHARGED BTZ BLACK HOLE (U = OO) 

A. The metric and the vector potential one-form 

For uj — oo, the corresponding three-dimensional charged Brans-Dicke theory is precisely the one provided by three- 
dimensional Einstein-Maxwell theory [2Cl I23L [25| . Then the general metric in Eq. J2]) , the <fi field in Eq. ((6]) , and the 
vector potential one-form ((SJ) reduce to the following 

*• . -(*-„- «f ,„ (f ) ) + g - — j*. _ _ + , (24, 

e" 20 = 1, (25) 

A = h Ql *{j) dT ' (26) 

where I is the AdS length, related to the cosmological constant by 2A 2 = l~ 2 , M is the mass, and Q is the electric 
charge. In this solution the scalar field 4> ls trivial. Next, in Fig. [TJ we show the Carter-Penrose diagram of the 
three-dimensional charged black hole. Despite there being a large and rich structure, the outside region in I is the 
relevant one, thus our foliation is from the bifurcation point to R — 00, and our boundary conditions reflect properties 
of the outer horizon at R = i?h, and of the spacelike infinity at R — 00. 



B. Canonical formalism 



Given Eqs. lf2"Tj) - ([23"ll . the three-dimensional action (U), now an Einstein-Maxwell action, becomes, excluding surface 
terms, 

S[A, R, T, A, R, f ; N, N r , $] = J dt J dr {-2N- x kR + 2N- 1 N r R'k + 2N- 1 (N r )'AR + 2N- 1 N r A'R 

^N^i^fA'R' + 2NA~ 2 A'R' - 2NA~ 1 R" + 4N\ 2 AR 
+ttN- 1 A- 1 r(t -*') 2 | ) (27) 

where 2A 2 = l~ 2 and I is the AdS length, which then yields the quantity a is a 2 = 2\ 2 = l~ 2 . In Eq. f27j) a ' means 
derivative with respect to time t and a ' means derivative with respect to r, and where all the explicit functional 
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FIG. 1: The Carter-Penrose diagram for the charged BTZ case, where is the outer horizon radius, Ri is the inner horizon 
radius, and the double line is the timelike singularity at R — 0. The static region I, from the bifurcation point to spacelike 
infinity R = oo, is the relevant region for our analysis. The regions II and IT are inside the outer horizon, beyond the foliation 
domain. Region I' is the symmetric of region I, towards left infinity. 



dependences are omitted. From the action 1|27|1 one determines the canonical momenta, conjugate to A, R, and T 
respectively, 

P A = -2N' 1 (it-R'N r \ , (28) 
P R = -2iV- 1 {A-(AiV r )'} , (29) 
Pr = 2?riV" 1 A" 1 i? (f - . (30) 



By performing a Legendre transformation, we obtain 

U = N \ —PrPa + — PrAiT 1 - 2A~ 2 A'P' + 2A~ 1 R" - 4A 2 Apl + N r {PrR 1 - P A A - TP r } + $ {-P r } 

^ 2 47T J 

= NH + N r H r + $G, (31) 

defining implicitly in this way the constraints H and H r , and G, and where the new Lagrange multiplier is <& = N r T. 
Here additional surface terms have been ignored, as for now we are interested in the bulk terms only. The action in 
Hamiltonian form is then 



S[A, R, T, P A , Pr, P r ; N, N r , #] = J dt dr {p A A + PrR + P r f - NH - N r H r - I>g} 



(32) 



The equations of motion are 

1 



A = --NP R + (N r AY, (33) 

R = -^NP A - N r R' , (34) 

Pr = 4X 2 NA- (2N , A- 1 Y + (N r PR,y + —P^NAR- 2 , (35) 

47T 

P A = 4\ 2 NR-2N'A~ 2 R' + N r (P A )' - —P^NR- 1 , (36) 

47T 

f = J-PrTVAp- 1 + (N r T)' + $' , (37) 
2n 

P r = iVPf. (38) 
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In order to have a well defined variational principle, we need to eliminate the surface terms of the original bulk action, 
which render the original action itself ill defined when one seeks a correct determination of the equations of motion 
through variational methods. These surface terms are eliminated through judicious choice of extra surface terms 
which should be added to the action. The action l132[) has the following extra surface terms, after variation 

Surface terms = (-N r P R 5R + N r ASP A + 2NA- 2 R'6A - 2NA- 1 6R / + 2N'A- 1 SR + N r T5P T + <t>5P r ) . (39) 

In order to evaluate this expression, we need to know the asymptotic conditions of each of the above functions 
individually, which are functions of (t,r). 
Starting with the limit r->0we assume 

(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 

Note that there are time dependences on the left hand side of the falloff conditions and that there are no such 
dependences on the right hand side, in the lower orders of the expansion in r. This apparent discrepancy stems from 
the fact that there is in fact no time dependence in the lower orders of the majority of the functions, but there may 
still exist such dependence for higher orders. Nevertheless, terms such as Rq are functions, independent of (t,r), thus 
constant, but undetermined. Their variation makes sense, as we may still vary between different values for these 
constant functions. With these conditions, we have for the surface terms at r = 



A(f, r) 


= A + O(r 2 ), 


R(t,r) 


= R Q + R 2 r 2 + 0(r 4 ) , 


PA{t,r) 


= 0(r 3 ), 


P R (t,r) 


= 0(r), 


N(t,r) 


= Ni(t)r + 0(r 3 ) , 


N r (t,r) 


= 0(r 3 ), 


r(t,r) 


= 0(r), 


Pr(t,r) 


= 2vrQ + 2irQ 2 r 2 + 0(r 4 ) , 


®(t,r) 


= $oW+0(r 2 ). 



Surface terms| r=0 = -2A r i A 1 SR - <M (2ttQ ) 



For r — » 00 we have 



A(t,r) 


= Ir- 1 + l 3 r](t)r- 3 + O x (r- 5 ) , 


R(t,r) 


= r + fptyr- 1 +0°°(r- 3 ), 


PA(t,r) 


= 0°°(r- 2 ), 


P R (t,r) 


= 0°°(r- 4 ), 


N(t,r) 


= i2(t,r) / A(t,r)- 1 (JV + (t)+0 00 (r- B )) ) 


N r (t,r) 


= 0°°(r- 2 ), 


T(t, r) 


= 0°°(r- 2 ), 


Pr(t,r) 


= 2iTQ + (t) + 0° a (r- 2 ), 


$(t,r) 


= ^ + (t)+0°°(r- 1 ). 



These conditions imply for the surface terms in the limit r — > 00 

Surface termsl^^ = 26{M+(t))N+ + $ +< 5 (2irQ+) , 
where M+(t) = 2(rj(t) + 2p{t)). So, the surface term added to (J32j is 

dt (^RoNiAo 1 - N+M+ + $0 (27tQ ) - $+ (2nQ+)^ 

What is left after varying this last surface term and adding it to the varied initial action (see Eq. (|32| ) is 

dt OiRodiN^o 1 ) - M+SN+ + 2ttQo5$>o - 2ttQ + 5$+ : 



SdT, 



A,R, Q ,Q+;N, $ ,$+ 



(49) 



(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 



(59) 



(60) 



(61) 
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We choose to fix -/ViA 1 and $o on the horizon (r = 0), and N+ and $o at infinity. These choices make the surface 
variation (|6lj) disappear. The term A^A^" 1 is the integrand of 

n a (h)n a (t 2 ) = -coHhQf dt N^Ao 1 (t)j , (62) 

which is the rate of the boost suffered by the future unit normal to the constant t hypersurfaces defined at the 
bifurcation circle, i.e., at r — » 0, due to the evolution of the constant t hypersurfaces. By fixing the integrand we are 
fixing the rate of the boost, which allows us to control the metric singularity when r — > Finally, we also fix <J> 
and <£>+ at r — and at infinity, respectively. 



C. Reconstruction, canonical transformation, and action 

In order to reconstruct the mass and the time from the canonical data, which amounts to making a canonical 
transformation, we have to rewrite the form of the solutions of Eqs. (J24j) - (J26j) . However, in addition to reconstructing 
the mass as done in [3] , we have to consider how to reconstruct the charge from the canonical data. We follow Kuchaf 
for this reconstruction. We concentrate our analysis on the right static region of the Carter-Penrose diagram. 

Developing the Killing time T as function of (t, r) in the expression for the vector potential (|26|) . and making use 
of the gauge freedom that allows us to write 

A = ^-Qln(?f)dT + dfi, (63) 



27r \l, 

where r) is an arbitrary continuous function of t and r, we can write the one-form potential as 

A = (i Qln (y) r '^')^ + (^ ln (f) T ^)*- ^ 

Equating the expression l|64p with Eq. (|23l> , and making use of the definition of Pr in Eq. (|30| , we arrive at 

Pr = 2ttQ. (65) 

In the static region, we define F as 

= ^-M-i-P r * In (!), (66) 

where we define the i?* as the (large) radius of a finite frontier, which will serve as a renormalization of the asymptotic 
properties of the future canonical coordinate M(t, r) (see subsection III Aft , This frontier will be used throughout the 
present section. We now make the following substitutions 

T = T(t,r), R = R(t,r), (67) 

into the solution (f24]) . getting 

ds 2 = -(FT 2 - F- x R 2 )dt 2 + 2(-FT'f + F- 1 R'R)dtdr + (-F(T'f + F" 1 R 2 ) dr 2 + R 2 d(p 2 . (68) 

This introduces the ADM foliation directly into the solutions. Comparing it with the ADM metric f2"Tj) . written in 
another form as 

ds 2 = —(N 2 - A 2 (N r ) 2 ) dt 2 + 2A 2 N r dtdr + A 2 dr 2 + R 2 dip 2 , (69) 
we can write a set of three equations 

A 2 = -F{T' f + F-\R') 2 , (70) 

A 2 iV r = —FT'f + F~ X R'R , (71) 

N 2 -A 2 (N r ) 2 = FT 2 — F~ X R 2 . (72) 
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The first two equations, Eqs. ((70) and Eq. fTTjl . give 



N r = -FTT + F-WR 

-F(T') 2 + F- 1 (R') 2 1 ' 

This one solution, together with Eq. (|70p. give 

R'f -T'R 

N = - . 74 

^-F{T'Y +F~ 1 (R>y 

One can show that N(t,r) is positive (see [13]). Next, putting Eqs. l(73|) - ((74) . into the definition of the conjugate 
momentum of the canonical coordinate A, given in Eq. lj28j) . one finds the spatial derivative of T(t, r) as a function 
of the canonical coordinates, i.e., 

— T' = -P^APa . (75) 

Later we will see that —T' = Pm, as it will be conjugate to a new canonical coordinate M. Following this procedure 
to the end, we may then find the form of the new coordinate M(t, r), as a function of t and r. First, we need to know 
the form of F as a function of the canonical pair A , R. For that, we replace back into Eq. ((70) the definition, in Eq. 
((75) . of T', giving 



Equating this form of F with Eq. l(66|) , we obtain 

£-£*'■'(£)-'■ (77) 

where F is given in Eq. Ip76|) , We thus have found the form of the new canonical coordinate, M. It is now a 
straightforward calculation to determine the Poisson bracket of this variable with Pm = —T' and see that they are 
conjugate, thus making Eq. fTB) the conjugate momentum of M, i.e., 

Pm = -P^APa • (78) 



The other natural transformation has been given in (J65J) , namely the one that relates Pr with Q through Pr = 2nQ. 
As M, Q is also a natural choice for canonical coordinate, being another physical parameter of the solution pi). The 
one coordinate which remains to be found is Pq, the conjugate momentum to the charge. It is also necessary to find 
out the other new canonical variable which commutes with M, Pm, and Q, and which guarantees, with its conjugate 
momentum, that the transformation from A, R, T, to M, Q, and the new variable is canonical. Immediately it 
is seen that R commutes with M, Pm, and Q. It is then a candidate. It remains to be seen whether Pr also 
commutes with M, Pm, and Q. As with R, it is straightforward to see that Pr does not commute with M and 
Pm, as these contain powers of R in their definitions, and {R(t,r), Pfl(i, r*)} = S(r — r*). So, rename the canonical 
variable R as R — R. We have then to find a new conjugate momentum to R which also commutes with M, Pm, 
and Q, making the transformation from {A, R, T; P\, Pr, Pr } — > {M, R, Q; Pm, Pr, Pq } a canonical one. The 
way to proceed is to look at the constraint H r , which is called in this formalism the super- momentum. This is the 
constraint which generates spatial diffeomorphisms in all variables. Its form, in the initial canonical coordinates, is 
H r = —AP A + PrR' — TPp. In this formulation, A is a spatial density, R is a spatial scalar, and V is also a spatial 
density. As the new variables, M , R, and Q, are spatial scalars, the generator of spatial diffeomorphisms is written 
as H r = PmM' + -PrR' + PqQ', regardless of the particular form of the canonical coordinate transformation. It is 
thus equating these two expressions of the super-momentum H r , with M, Pm, and Q written as functions of A, R, T 
and their respective momenta, that gives us one equation for the new Pj^ and Pq. This means that we have two 
unknowns, Pj^ and Pq, for one equation only. This sugests that we should make the coefficients of R' and equal 
to zero independently. This results in 

P R = P R - 2\ 2 F- 1 AP A R + F -1 A _1 P A ii" - P _1 A" 2 P A A'P' 

-F^A^P'^R' + —P$R- 1 F- 1 AP A , (79) 

P Q = -27rr+ip r ln('-^')p- 1 APA. (80) 
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It can be shown that Pq commutes with the new P R and with the rest of the new coordinates, except with Q. We 
have now all the canonical variables of the new set determined. For completeness and future use, we write the inverse 
transformation for A and Pa, 



A = ({R'fF- 1 - P 2 M F) 
Pa = 2FP M ((R"fF- 1 - 
In summary, the canonical transformations are the following, 

R = P, 



PltF) 



(81) 
(82) 



R 



1 



M = — PA In — \-F 



Aw 



R 



R* 



Q = 
P R = 



r 2 

Pr 
2ir ' 

P R - 2X 2 F- 1 AP A R + F- x A~ l P A R" 
1 



P- 1 A" 2 P A A'P / 



-F^A^P^R' 



1 



P M = -P AP y 
Po -- 



8tt 



R 



PZR-^-F^AP, 



A - 



2tiT + -P r In — P _1 AP 



A • 



(83) 



It remains to be seen that this set of transformations is in fact canonical. In order to prove that the set of equalities 
in expression 1(83]) is canonical we start with the equality 



P A 6A + P R 5R + P r ST - P M 5M - P R 6R - P Q c 



SR\n 

S 



2R' + AP A 



2R 1 - AP A 
rp r + AP A + 



R' In 



2R' - 


-AP A 


) 


2R' - 


f AP A 





(84) 



We now integrate expression (|84|) in r, in the interval from r = to r — oo. The first term on the right hand side of 
Eq. dm vanishes due to the falloff conditions (see Eqs. fl0 |l -(|i8 |i and Eqs. ([50j) - ((SHJ) ) - We then obtain the following 
expression 



/>oo />oo 

/ dr (Pa 5 A + P B SR + P r ST) dr (P M SM + P R SR + P Q SQ) = Suj [A, R, T, P A , P r ] , 

Jo Jo 



(85) 



where Suj [A, R, T, Pa, Pr] is a well defined functional, which is also an exact form. This equality shows that the 
difference between the Liouville form of {P, A, T; Pr, Pa, Pr} and the Liouville form of {R, M, Q; P r , Pm, Pq} is 
an exact form, which implies that the transformation of variables given by the set of equations f83|) is canonical. 

Armed with the certainty of the canonicity of the new variables, we can write the asymptotic form of the canonical 
variables and of the metric function F(t, r). These are, for r — > 

(86) 
(87) 



F(t,r) 


= ARlA^r 2 + 


0(r 4 ), 


R(t,r) 


= R + R 2 r 2 - 


H 0(r 4 ) , 


M (t, r) 


= M Q + M 2 r 2 


+ 0(r 4 ), 


Q{t,r) 


= Qo + Q2r 2 - 


fO(r 4 ), 




= 0(r), 




P M (t,r) 


= 0(r), 




PQ(t,r) 


= 0(r), 





with 



= l-'Rl--KQlhx[ ^ 



P, t 



M 

Mi = 2r z R Q R 2 - 27iTn I ^ ) Q Q 2 - 



Rr 



)lR 2 R 1 - 4A 2 R\ 



(89) 
(90) 
(91) 
(92) 



(93) 
(94) 
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For r — > oo, we have 



F(t,r) 


= 2X 2 r 2 ~2{r 1 (t) + 2p(t)) + O ao (r- 2 ), 


(95) 


R(t,r) 


= r + (2X 2 )- 1 p(t)r- 1 +0°°(r- 3 ), 


(96) 


M(t,r) 


= M + (t)+0°°(r- 2 ), 


(97) 


Q(t,r) 


= Q + (t) + 0™(r- 2 ), 


(98) 




= 0°°(r- 4 ), 


(99) 


P M (t,r) 


= 0°°(r- 5 ), 


(100) 


PQ(t,r) 


= O 00 ^ 2 ), 


(101) 



where M + (t) = 2(rj(t) + 2p(t)), as seen before in the surface terms (see Eq. l(60|)). Note here that without the 
treatment of the divergence in subsection III AI the function M(t,r) would diverge. 

We are now almost ready to write the action with the new canonical variables. It is now necessary to determine 
the new Lagrange multipliers. In order to write the new constraints with the new Lagrange multipliers, we can use 
the identity given by the space derivative of M, 



M' = -A - 



R'H + ip A (H r - TG) - -^PrGln ( £ ) A 

Z Z7T V R 



(102) 



Solving for H and making use of the inverse transformations of A and Pa, in Eqs. lf8lj) and ([82]) . we get 

M'F-iR' + FP M P R + 2ttQQ'R'F- 1 In 



H 



(F-i(R') 2 -^) 



H r = P M M' + P R R' + P Q 
G = -2txQ' . 



(103) 

(104) 
(105) 



Following Kuchaf [10], the new set of constraints, totally equivalent to the old set H(t,r) = 0, H r (t,r) = 0, and 
G = 0, outside the horizon points, is M'(t,r) = 0, P R (t,r) = 0, and Q'(t,r) = 0. By continuity, this also applies on 
the horizon, where F(t, r) — 0. So we can say that the equivalence is valid everywhere. The new Hamiltonian, which 
is the total sum of the constraints, can now be written as 



NH + N r H r + $G = N M M' + N R Pt> + N Q Q' 



R 



(106) 



In order to determine the new Lagrange multipliers, one has to write the left hand side of the previous equation, Eq. 
(|106p . and replace the constraints on that side by their expressions as functions of the new canonical coordinates, 
spelt out in Eqs. (| 103[) - (j 105p . After manipulation, one gets 



N 



M 



NF~ 1 FL' 



(F- l {R'f -FPf 



- + N r P M , 



M! 



R 



N 



NFP M 



{F-^(R'f-FP 2 M y 
2TrNQR'F- 1 \n(^- 



- + N r R' , 



(F-1(K 



'\2 



FPli) 



N t Pq - 2tt$ . 



(107) 
(108) 

(109) 



Using the inverse transformations Eqs. ([8l"Tl - lj82"j) . and the identity R = R, with Pr — 2ttQ, we can write the new 
multipliers as functions of the old variables 



N M 

N R 



-NF^R'A- 1 + -N'F^KPa 

1np a + n t r', 



= NPrR'F^A- 1 In ( ^- ) - 27riVT + ^N^rF- 1 In ( -^-] AP A - 2tt$ , 



R* 



R* 



(110) 
(111) 
(112) 



12 



allowing us to determine its asymptotic conditions from the original conditions given above. These transformations 
are non-singular for r > 0. As before, for r — > 0, 

N M (t,r) = -l Nl (t)A R^ +0(r 2 ), (113) 

N R (t,r) = 0(r 4 ), (114) 

N Q (t,r) = ii? 2 - 1 ln^)g A 7Vi(t)-2^ (i)+O(r 2 ), (115) 

and for r-»oowe have 



N M (t,r) = -N + (t)+0°°(r- 2 ), (116) 
N R {t,r) = O 00 ^- 1 ), (117) 
N Q (t,r) = -27rl- + (t) + 00 (r- 1 ). (118) 

The conditions (| 1 13jl - (| 1 18|) show that the transformations in Eqs. f| 1 10(1 - (| 1 1 1|) are satisfactory in the case of r — ► oo, 
but not for r — > 0. This is due to fact that in order to fix the Lagrange multipliers for r — ► oo, as we are free to do, 
we fix N + (t), which we already do when adding the surface term 

dtN+M+ (119) 

to the action, in order to obtain the equations of motion in the bulk, without surface terms. The same is true for $ + . 
However, at r = 0, we see that fixing the multiplier N M to values independent of the canonical variables is not the 
same as fixing A^Aq 1 to values independent of the canonical variables. The same is true of the fixation of with 
respect to $o- We need to rewrite the multipliers N M and for the asymptotic regime r — > without affecting 
their behavior for r — ► oo. In order to proceed we have to make one assumption, which is that the expression given 

in the asymptotic condition of M(t,r), as r — * 0, for the term of order zero, Mo = 1~ 2 Rq — 7rQg In ("r 1 )) defines 

Ro as a function of Mo and Qo, and Rq is the horizon radius function, Ro = R^(Mq,Qq). Also, we assume that 
M > M crit (Q ), where 



M ctit (Q ) = ^ ^1 - In ^ j j . (120) 

With these assumptions, we are working in the domain of the classical solutions. We can immediately obtain that the 
variation of Rq is given in relation to the variations of M and Qo as 



SRo = {2l- 2 R - irQlRv 1 ) UMq + 2^Q In SQo J , (121) 

This expression will be used when we derive the equations of motion from the new action. We now define the new 
multipliers N M and as 



N M = -N M 



(l-g) + 2g (2r 2 R a - t:Q 2 R^) , (122) 



NQ = N M 2gQ ln(^J (2l~ 2 R — ttQqRq 1 ) 1 — , (123) 

where g(r) = 1 + 0(r 2 ) for r — > and g(r) = 0°°(r~ 5 ) for r — ► oo. The new multipliers, functions of the old 
multipliers N M and , have as their properties for r — > 

A> M (t,r) = A> M (i) + O(r 2 ), (124) 
A Q (t,r) = 2vrl>o(t) + 0(r 2 ), (125) 



and as their properties for 



N M {t,r) = A>+(i) + 0°°(r- 2 ), (126) 
N Q (t,r) = 2 7 rl> + (t) + oo (r- 1 ). (127) 
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When the constraints M' = = Q' hold, Nq is given by 



N, 



M 



(128) 



With this new constraint N M , fixing A^Aq 1 at r = or fixing Nq 1 is equivalent, there being no problems with N R , 
which is left as determined in Eq. (| 108|1 . With respect to the same happens, i. e., fixing the zero order term of 
the expansion of for r -> is the same as fixing $o ■ At infinity there were no initial problems with the definitions 
of both N M and N Q . 

The new action is now written as the sum of Ss, the bulk action, and San, the surface action, 



S 



M, R, Q, Pm,P R , Pq; N M , N R , N q 



dt J dr [P M M + P R R 
N Q Q' - N R P R + N M (1 - g) M' 
N M 2g (2r 2 R Q - nQlRo 1 

dt | (2R N^ - N+M+ 



M' - Q'Q In 



Rq 



2tt ($ Qo 



} • (129) 



The new equations of motion are now 



M 
R 

Q 

Pm 

p r 

Pq 
m' 

p r 

Q' 



o, 

iV R , 

o, 

(N M y, 
o, 

(N Q y, 

o, 

0, 
0. 



(130) 
(131) 
(132) 
(133) 
(134) 
(135) 
(136) 
(137) 
(138) 



where we understood N M to be a function of the new constraint, defined through Eq. (|122|) and as the function 
of the new constraint defined through Eq. (|123j) . The resulting boundary terms of the variation of this new action, 
Eq. (|129|) . are, first, terms proportional to SM, 5R, and SQ on the initial and final hypersurfaces, and, second, 



dt ( 2R 6N^ - M+SN^ 



2tt [Qo5§o - Q+SQi 



(139) 



To arrive at (|139|1 we have used the expression in Eq. (|121|) . The action in Eq. (|129|) yields the equations of motion, 
Eqs. I|130p - I|138p . provided that we fix the initial and final values of the new canonical variables and that we also fix 
the values of Nq 1 and of N+, and of $o and $+. Thanks to the redefinition of the Lagrange multiplier, from N M to 
N M , the fixation of those quantities, N^f and N + , has the same meaning it had before the canonical transformations 
and the redefinition of iV M . The same happens with $ and $ + . This keeping of meaning is guaranteed through 
the use of our gauge freedom to choose the multipliers, and at the same time not fixing the boundary variations 
independently of the choice of Lagrange multipliers, which in turn allow us to have a well defined variational principle 
for the action. 



D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion 
(|130p - l|138p allow us to write M and Q as independent functions of space, r, 



M(t,r) = m(t), 
Q(t,r) = q(t). 



(140) 
(141) 
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The reduced action, with the constraints taken into account, is then 



m,p m ,q, p q 



N M 



dt p m m + p q q h , 



where 



Pm 



Pq 



drP, 



M ■ 



drP c 



and the reduced Hamiltonian, h, is now written as 



h(m, q; t) = -2R h N ( f + N+m + 2nq ( 



$ 



(142) 

(143) 
(144) 

(145) 



with i?^ being the horizon radius. We also have that m > M crit (q), according to the assumptions made in the 
previous subsection. Thanks to the functions Nq 1 ^), N+(t), $o(t), and $+(t) the Hamiltonian h is an explicitly time 
dependent function. The variational principle associated with the reduced action, Eq. (jT42j) , will fix the values of 
m and q on the initial and final hypersurfaces, or in the spirit of the classical analytical mechanics, the Hamiltonian 
principle fixes the initial and final values of the canonical coordinates. The equations of motion are 



q 

Pm 
Pq 



2n {2 qln (jj^j (2l~ 2 R h - TrQgi^ 1 )" 1 ^ + $0 - $+} 



(146) 
(147) 
(148) 

(149) 



The equation of motion for m, Eq. (| 146|1 . is understood as saying that m is, on a classical solution, equal to the mass 
parameter M of the solution, Eq. ([24|) . The same goes for the function q, where Eq. I|147p implies that q is equal 
to the charge parameter Q on a classical solution, Eq. (|24| . In order to interpret the other equation of motion, Eq. 
(|148jl . we have to recall that from Eq. f78|) one has Pm = —T', where T is the Killing time. This, together with the 
definition of p m , given in Eq. (| 143(1 . yields 



Pm 



Tn-T, 



(150) 



where To is the value of the Killing time at the left end of the hypersurface of a certain t, and T + is the Killing time 
at spatial infinity, the right end of the same hypersurface of t. As the hypersurface evolves in the spacetime of the 
black hole solution, the right hand side of Eq. <|147|) is equal to T — T + . Finally, after the definition 



Pq = £0 £4 



(151) 



obtained from Eqs. ([64]) . |80|) . and <|144[1 . Eq. I|149p gives in the right hand side £0 — £+, which is the diference of the 
time derivatives of the electromagnetic gauge £(t, r) at r = and at infinity. 



E. Quantum theory and partition function 



The next step is to quantize the reduced Hamiltonian theory, by building the time evolution operator quantum 
mechanically and then obtaining a partition function through the analytic continuation of the same operator [l|- 
0| . The variables m and q are regarded here as configuration variables. These variables satisfy the inequality 
m > M^t (q). The wave functions will be of the form ip(m, q), with the inner product given by 

(V , ,X)= / jueZmdq^x, (152) 
J A 

where A is the domain of integration defined by m > M crit (q) and /u(m, q) is a smooth and positive weight factor for 
the integration measure. It is assumed that fi is a slow varying function, otherwise arbitrary. We are thus working in 
the Hilbert space defined as J$? := L 2 (A; [idmdq). 
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The Hamiltonian operator, written as h(t), acts through pointwise multiplication by the function h(m, q;i), which 
on a function of our working Hilbert space reads 



h(t)tp(m, q) = h(m, q; t)ip(m, q) . 



(153) 



This Hamiltonian operator is an unbounded essentially self-adjoint operator. The corresponding time evolution 
operator in the same Hilbert space, which is unitary due to the fact that the Hamiltonian operator is self-adjoint, is 



exp 



I2 



dt' h(i') 



ti 



This operator acts also by pointwise multiplication in the Hilbert space. We now define 

fl-2 



T 



6 := 



-1) 



tl 



dtN+(t) , 



dtN™(t), 



dt$ Q (t), 



dt$+(t) . 



(154) 

(155) 
(156) 
(157) 
(158) 



Using l|145p . |154p - (|158p we write the function K, which is in fact the action of the operator in the Hilbert space, as 

K(m;T,Q,E ,E + ) = exp [-imT + 2 i i? h 9 - 2tt iq(E+ - Ho)] . (159) 

This expression indicates that Kfeiti) depends on t\ and t 2 only through the functions T, 6, H , and 3+. Thus, 
the operator corresponding to the function K can now be written as K(T, 9, H , H + ). The composition law in time 
K(ta;t2)K(t2',ti) — K(t 3 ;ti) can be regarded as a sum of the parameters T, 6, Ho, and H + inside the operator 
K(T, 0, H , H + ). These parameters are evolutions parameters defined by the boundary conditions, i.e., T is the 
Killing time elapsed at right spatial infinity and 8 is the boost parameter elapsed at the bifurcation circle; Ho and 
H + are line integrals along timelike curves of constant r, and constant angular variables, at r = and at infinity. 



F. Thermodynamics 

We can now build the partition function for this system. The path to follow is to continue the operator to imaginary 
time and take the trace over a complete orthogonal basis. Our classical thermodynamic situation consists of a three- 
dimensional spherically symmetric charged black hole, asymptotically anti-de Sitter, in thermal equilibrium with a 
bath of Hawking radiation. Ignoring back reaction from the radiation, the geometry is described by the solutions in 
Eqs. l(24|) - (|26l) . Thus, we consider a thermodynamic ensemble in which the temperature, or more appropriately here, 
the inverse temperature j3, and the electric potential (j> are fixed. This characterizes a grand canonical ensemble, and 
the partition function Z(f3, <j>) arises naturally in such an ensemble. To analytically continue the Lorentzian solution 
we put T = —i/3, and 9 — 2tti, this latter choice based on the regularity of the classical Euclidean solution. We also 
choose Ho = and H + = if3(f>. 

We arrive then at the following expression for the partition function 

Z(J3, $) = Tr \k{-%P, -2m, 0, ify) . (160) 

From Eq. (| 159|) this is realized as 

Z{P4) = / fidmdq exp \-(3(m - 27rq0) + 47ri? h ] (m|m) . (161) 

J A 

Since (m|m) is equal to 6(0), one has to regularize Eq. (|16ip . Following the procedure developed in the Louko- Whiting 
approach Q, this means regularizing and normalizing the operator K beforehand. This leads to 

2xen(/3, $) = Af J (i dmdq exp [-/3(m - 2nq$) + 47ri? h ] , (162) 
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where A/" is a normalization factor and A is the domain of integration. Provided the weight factor /i is slowly varying 
compared to the exponential in Eq. I|162p , and using the fact that the horizon radius is a function of m and q, 
the integral in Eq. (|162|) is convergent. Changing integration variables, from m to i? n , keeping q, where 



2 D 2 



7rq In 



the integral Eq. I|162p becomes 



-2ren(/3, 4>)=N I ndR-^dq exp(-Q , 



(163) 



(164) 



where A' is new the domain of integration after changing variables, and the function 2*(i? n ), a kind of an effective 
action (see [IB]), is written as 



I,(R h ) :=P r^R h -7rq 2 ln 



R* 



(165) 



The new domain of integration, A' , is defined by the inequalities < R^ and q 2 < 2i? 2 j 7r l l 2 . The new weight 
factor Jl includes the Jacobian of the transformation. Since the weight factor Jl is slowly varying, we can estimate the 
integral of Z re n(/3, </>) by the saddle point approximation. For that we have to calculate the critical points. Firstly 
one finds the value of q for which 



dI*(R h ,q) 
<9q 



The value is 



In 



= 0. 



Replacing this value in l|165p . we obtain 

I,(R h ,q*(R h )) = (3[r 2 R\ 



AirRi- 



Deriving expression (| 168|) with respect to R^ and making it zero, we obtain 

n -2 



gUgfa) 

dR h 



= 2f3r 2 R h - 0<j> 2 



In 



^h 



4tt = 0. 



which implies 



2f3r 2 R? - ffi 



In 



^h 



AwR-h = 0. 



(166) 



(167) 



(168) 



(169) 



(170) 



The solutions R^ to the equation (| 1 T0[) gives us the pair (R^, q*) of critical points, with q* given in Eq. (|167|1 . where 
= i?^. Of these, the classical solution comes from the one critical pair that is also the global minimum of the 
effective action /*, denoted by (i?^,q + ), depending on the interval of values to which (f) belongs. The renormalized 
partition function is then 



2ren(/5» = Pexp 



-0 (r 2 (Rp 2 -n( q +) 2 \n(^ 



27rq 4 



4vri?+ 



(171) 



where P is a slowly varying prefactor and {Rt, q + ) is the global minimum of the effective action (j 165[) . In the domain 
of integration the dominating contribution comes from the vicinity of R^ = R^. We now write the logarithm of 2 Ten 



ln(Z ren ) = lrxP - 01 r 2 (Rp* - Tr(q + y In 



5 

JL 



4 ^h' 



(172) 
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By ignoring the prefactor's logarithm, which closer to R^ is less relevant, we are able to determine the value of m 
at the critical point, where we find that it corresponds to the value of the mass of the classical solutions of the black 
holes (see Eq. (|24|) ). Thus, when the critical point dominates the partition function, we have that the mean energy 
(E) is given by 

(E) = -^lnZ re n « 2A 2 (i?+) 2 - 7r(q+) 2 ln (^j = m+ , (173) 
where m + is obtained from Eq. l|163j) evaluated at R^ . The thermal expectation value of the charge is 

(G> =/3" 1 4 lnZ ren~2 7 rq+. (174) 
The temperature of the black hole, T = /3 _1 , is 

T = i- (4A 2 (i?+) - 7r(q+) 2 (i?+)- 1 ) , (175) 

where q + is the function in 1|167|1 evaluated at Rt. If the maximum value of the charge q + is chosen, i.e., q 2 — 
2R?tt~ 1 1~ 2 , which is the value of the charge of the extreme solution of the black hole, and is replaced into (j 1 T5[> . then 
the temperature is null, as expected from an extreme solution. It can be shown that 9m + /d/3 < 0, which through the 
constant 4> heat capacity — —f3 2 (d (E) /d(3) tells us that the system is thermodynamically stable. The entropy is 
given by 

S=(l- p (In Z ren ) « 47ri?+ . (176) 

This is the entropy of the charged BTZ black hole [lllHS] (see also (HI)- This entropy includes the extreme solution, 
as the specific value of the charge is irrelevant for the determination of the entropy, provided we are in the domain of 
validity of the approximation and q 2 < 2R 2 ^ir~ 1 l~ 2 . 



IV. HAMILTONIAN THERMODYNAMICS OF THE GENERAL RELATIVISTIC CYLINDRICAL 
DIMENSIONALLY REDUCED CHARGED BLACK HOLE (U) = 0) 

A. The metric, the scalar 0, and the vector potential one-form 



For u> = 0, the corresponding three-dimensional charged Brans-Dicke theory is obtained from the cylindrical di- 
mensionally reduced black hole of four-dimensional general relativity, with a Maxwell term (2lL I22I |25| . Then general 
metric in Eq. the <f> field in Eq. iJH), and the vector potential one-form ([8]) reduce to the following 



ds 2 = 
e- % * = aR 



(ai?) 2 



M 



Q 2 



2(aR) 16(ai?) 5 



dT 2 



dR 2 



{aRf 



M 



+ R 2 d<p 2 , 



2{aR) ' 16(aR) 2 



A = 



(177) 

(178) 
(179) 



with M — 26, a = y ||A| = I l , and a — 4a, where the I is the AdS length. The Schwarzschild coordinates are 

again {T,R,ip}. Unlike the charged BTZ solution, this solution (|177p - l|179p has a metric function whose mass term 
depends on R. This behavior is similar to the Reissner-Nordstrom black hole metric function jll|. In Fig. [2] we show 
the Carter-Penrose diagram of the black hole solution for w = 0, where again i?^ is the horizon radius, given by the 
larger positive real root of 



{aR h f- l -M{aR h )+(Q-\ =0. 



(180) 



R = is the radius of the timelike curvature singularity, and R — oo is the spatial infinity. 
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FIG. 2: The Carter-Penrose diagram for the charged u) = case, where is the outer horizon radius, -Ri is the inner horizon 
radius, and the double line is the timelike singularity. 



B. Canonical formalism 



The action with lu = becomes, following Eqs. i[2l|) - ([23| . and excluding surface terms, 

S[A, R, T, A, R, f ; N, N r , $] = J dt dra | X 2 ANR 2 - N^RRA - ^N~ 1 AR 2 + N- 1 RR{AN r )' 

+N^N r RR'A - N- 1 N r ARR / (N r Y - N' 1 (N r ) 2 RR' A' + N~ 1 N r ARR' 
_lAr-i(7v' ) 2 A(i?') 2 - (A-^'RR'N - ^A'^R'^N - A^Rtf'N 

+ ~aN- 1 A~ 1 R 2 (t - $') 2 | , (181) 

where ' means derivative with respect to time t and ' is the derivative with respect to r, and where all the explicit 
functional dependences are omitted. Depending on the situation we use four different letters containing the same 
information but with slightly different numerical values. Thus, a, A, I and a are related by a 2 = |A 2 = Z~ 2 , and 
a — 4a, where I is the AdS length. From this action one determines the canonical momenta, conjugate to A, R, and 
r respectively 

Pa = -aN^R^R- R'N r } , (182) 
P R = -aN- 1 {p[A - (AN r )'\ + A[R - N r R'}} , (183) 
P r = aN~ 1 A~ 1 R 2 (f - $') . (184) 



By performing a Legendre transformation, we obtain 

H = N i—a~^R~ 1 P\Pn + ^a^AR^iPl + P$) + aA^RR" - aA- 2 RR'A' + ^A^P') 2 - A 2 AP 2 

+N r {P R R' - P\A - Pf r} + $ {-Pi } = NH + N r H r + $G , (185) 
with the new Lagrange multiplier being <f> := <I> — APT. The action in Hamiltonian form is then 

S[A, P, T, P A , P R , P r ; N, N r , $] = J dtj^ dr |p A A + P R R + P r f - NH - N r H r - $g| . (186) 
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The equations of motion are 

A 
R 

r 

Pn 



-Na^R-iPa + Na~ 1 AR~ 2 P^ + N r A , 



Na^P^R- 1 + N r R' , 
^NAR-ZPr + (N r T)'& , 

Na^P K P R R- 2 + Na^kiPl + F^iT 3 - ((A^/A^i?)' - (NajiA^R')' 
2Na\ 2 AR + (N r P R )' , 



P A = -^Na~ 1 Rr 2 {Pl + P^) - (NaYRR'A- 2 - ^Na(R') 2 A- 2 
P r = N r P^. 



NaX 2 R 2 



N r P A , 



(187) 
(188) 
(189) 

(190) 
(191) 
(192) 



In order to have a well defined variational principle, we need to eliminate the surface terms of the original bulk action, 
which render the original action itself ill defined for a correct determination of the equations of motion through 
variational methods. Through the choice of added surface terms one can achieve this elimination. The action (|186p 
has the following extra surface terms, after variation 



Surface terms = {a (-NA^RSR! + N' A" 1 RSR - N r P R 5R + N r ASP A + NRR' A~ 2 8A) 



+N r T6Pr + <S>SP r 



(193) 



In order to evaluate this expression, we need to know the asymptotic conditions of each of the functions of (t, r). 
Starting with the limit r^Owe assume 



A(t, r) 


= A (t)4 


-0(r 2 ), 




R(t, r) 


= flo(*H 


-R 2 (t)r 2 ^ 


-0(r 4 ), 


r(t,r) 


= 0(r), 






P\(t,r) 


= 0(r 3 ), 






Pn(jt,r) 


= 0(r), 






Pr(t,r) 


= Qo(*H 


- Q2{t)r 2 4 


-0(r 4 ), 


N(t,r) 


= Ni{t)r 


+ 0(r 3 ), 




N r {t 7 r) 


= Nl{t)r 


+ 0(r 3 ), 




$(t,r) 


= $o(*)4 


0(r 2 ). 





With these conditions, we have for the surface terms at r = 0, 



Surface terms 



r=0 



-ocNiRqA^SRo - 3>o5Qo 



In the same way, for r — > oo, 



A(t,r) 


= ?r- 1 + Z 3 ?7(t)r- 4 + oo (r- 5 ), 


R(t,r) 


= r + l 2 p(t)r- 2 + 0°°(r- 3 ), 


r(i, r) 


= O 00 ^ 2 ), 


P\(t t r) 


= 0°°(r- 2 ), 


Pn(t,r) 


= 0°°(r- 4 ), 


Pr(t,r) 


- Q+W + O 00 ^ 1 ), 


N(t,r) 


= i2(t,r) / A(t,r)- 1 (JV- + (t)+0 00 (r- B )) ) 


N r {t,r) 


= 0°°(r- 2 ), 




= +O oc '(r- 1 ). 



These conditions imply for the surface terms in the limit r — > oo, 



Surface terms^^ = a<5(A/+)Ar+ + &+5Q+ 



where M+(t) — ot(r){t) + 3p(t)). So, the surface term to be added to (|186|1 is 



Sdi 



A,R,Q ,Q + ;N,$ ,$ + = j dt ( -aR^NtA^ 1 - N+M+ + $> Q - $+Q+ 



(194) 
(195) 
(196) 
(197) 
(198) 
(199) 
(200) 
(201) 
(202) 

(203) 

(204) 
(205) 
(206) 
(207) 
(208) 
(209) 
(210) 
(211) 
(212) 

(213) 



(214) 
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What is left after varying this last surface term and adding it to the varied initial action (see Eq. (|186p ) is 

dt (^aB^SiNxAo 1 ) - M+SN+ + Q 5$ - Q+SQ+^J . (215) 

We choose to fix NiAq 1 and $o on the horizon, and N + and at infinity, which makes the surface variation (|2 1 5j) 
disappear. 

C. Reconstruction, canonical transformation, and action 

In order to reconstruct the mass and the time from the canonical data, which amounts to making a canonical 
transformation, we have to rewrite the general form of the solutions in Eqs. (|177|1 - (|179[1 . However, in addition to 
reconstructing the mass as done in [l[ , we have to consider how to reconstruct the charge from the canonical data. We 
follow Kuchaf [ioj for this reconstruction. We concentrate our analysis on the right static region of the Carter-Penrose 
diagram. 

Developing the Killing time T as function of (t, r) in the expression for the vector potential <|179[) . and making use 
of the gauge freedom that allows us to write 

A = oT x %dT + d£, (216) 
R 

where £(t, r) is an arbitrary continuous function of t and r, we can write the one-form potential as 

A = (a-^T' + e^dr+L-^T + ijdt. (217) 

Equating the expression l|217j) with Eq. (|23| , and making use of the definition of Pr in Eq. (| 184ft , we arrive at 

P r = Q. (218) 

In the right static region we define F as 

M P 2 

F(R(t, r)) = (aR(t, r)f - + ■ (219) 

Note that here, i?*, the large radius of a finite frontier, which can be used to renormalize the asymptotic properties of 
the future canonical coordinate M(t,r) (see subsections III Al and IIII C|l , is not needed; in fact bringing to the to = 
solution introduces complications at the r = frontier. Thus, in this section we do not mention i?*. Then, making 
the following substitutions 

T = T(t,r), R = R(t,r), (220) 

in the solutions above, Eqs. (jT77j) - (jT7Sj) , one has 

ds 2 = -{FT 2 - F- x R 2 )dt 2 + 2(-FT'f + F~ 1 R'R)dtdr + (~F(T'f + F" 1 R 2 ) dr 2 + R 2 d(p 2 . (221) 

This introduces the ADM foliation directly into the solutions. Comparing it with the ADM metric l(2"Tj) . written in 
another form as 

ds 2 = -(N 2 - A 2 (N r ) 2 ) dt 2 + 2A 2 N r dtdr + A 2 dr 2 + R 2 dip 2 , (222) 
we can write a set of three equations 

A 2 = -F(T') 2 + F-\R') 2 , (223) 

A 2 iV r = -FT'f + F^R'R, (224) 

N 2 -A 2 {N r ) 2 = FT^-F-^R 2 . (225) 

The first two equations, Eq. ([223]) and Eq. ([224| . give 

Nr = -F(T>) 2 + F-W (226) 
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This one solution, together with Eq. (|223[> . give 



R T — T R , „, 

N = , . 227 
^-F(T') 2 +F- l {R'f 

One can show that N(t, r) is positive. Next, putting Eq. <|226[1 and Eq. I|227p into the definition of the conjugate 
momentum of the canonical coordinate A, given in Eq. (| 182[) . one finds the spatial derivative of T(t, r) as a function 
of the canonical coordinates, i.e., 

- T' = oT 1 RT X P" 1 AP A . (228) 

Later we will see that — T' = Pm, as it will be conjugate to a new canonical coordinate M. Following this procedure 
to the end, we may then find the form of the new coordinate M(t, r), also as a function of t and r. First, we need to 
know the form of F as a function of the canonical pair A , R. For that, we replace back into Eq. I|223p the definition 
of X" , giving 

R'Y /P\ N 1 



F =( A -U|) ' (229) 



Equating this form of F with Eq. l|219j) , we obtain 

1 fa 2 2 a 2 R 2 \ , x 

where F is given in Eq. ((229]). We thus have found the form of the new canonical coordinate, M. It is now a 
straightforward calculation to determine the Poisson bracket of this variable with Pm = —T' and see that they are 
conjugate, thus making Eq. (|228|) the conjugate momentum of M, i.e., 

P M = a _1 i2~ 1 J , - 1 AP A . (231) 

The other natural transformation has been given, namely the one that relates Pr with Q through Pr = Q. As 
M, Q is also a natural choice for canonical coordinate, being another physical parameter of the solution l[2"4"|). The 
one coordinate which remains to be found is Pq, the conjugate momentum to the charge. It is also necessary to find 
out the other new canonical variable which commutes with M, Pm, and Q, and which guarantees, with its conjugate 
momentum, that the transformation from A, R, T, to M, Q, and the new variable is canonical. Immediately is 
it seen that R commutes with M, Pm, and Q. It is then a candidate. It remains to be seen whether Pr also 
commutes with M, Pm, and Q. As with R, it is straightforward to see that Pr does not commute with M and 
Pm, as these contain powers of R in their definitions, and {R(t,r), Pn(t,r*)} = S(r — r*). So rename the canonical 
variable R as R — R. We have then to find a new conjugate momentum to R which also commutes with M, Pm, 
and Q, making the transformation from {A, R, T; P\, Pr, Pr } — > {M, R, Q; Pm, Pr, Pq } a canonical one. The 
way to proceed is to look at the constraint H r , which is called in this formalism the super- momentum. This is the 
constraint which generates spatial diffeomorphisms in all variables. Its form, in the initial canonical coordinates, is 
H r = —AP A + PrR' — TPp. In this formulation, A is a spatial density, R is a spatial scalar, and T is also a spatial 
density. As the new variables, M, R, and Q, are spatial scalars, the generator of spatial diffeomorphisms is written 
as H r = PmM' + PrR' + PqQ', regardless of the particular form of the canonical coordinate transformation. It is 
thus equating these two expressions of the super-momentum H r , with M, Pm, and Q written as functions of A, R, T 
and their respective momenta, that gives us one equation for the new P^ and Pq. This means that we have two 
unknowns, sc. Pr and Pq, for one equation only. This sugests that we should make the coefficients of R' and Pf 
equal to zero independently. This results in 

P R = p R ~ 1& F ~ 1APaR ~ \ R ~ lAP * + F^PaRJ'A- 1 - F^A^PaA'BJ + (R'^F^A^PaR- 1 

-F^A^PaR' + -a~ 2 i?~ 3 J^F -1 APA , (232) 
P Q = -T-a~ 2 R- 2 F- l AP A Pr. (233) 

We have now all the canonical variables of the new set determined. For completeness and future use, we write the 
inverse transformation for A and Pa, 

A = ((R') 2 P _1 -PmFY , (234) 
Pa = oRFPm {{ti fF- 1 - P 2 M Fy k . (235) 
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In summary, the canonical transformations are 
R = R, 



M = 



R 



1 

2 
Pr 



Pi 



3a' 
"32 



F^APaR - -R^APa + F -1 PaE"A _1 - F _1 A- 2 P A A'P' + (P') 2 P" 1 A" 1 PaP _1 



-F^A^PiR' 



1 



-■2r-3 P 2 f -i ap 



Pm 

Pq 



a~ 1 R~ 1 F- L AP A , 

-r - a" 2 P _2 F _1 AP A Pr . 



(236) 



In order to prove that the set of equalities in expression (|236|1 is canonical we start with the equality 



P A <5A + Pi?<SP + P r 5r - PmSM - P R 8R - P Q SQ 



-aR5R\n 
2 



aPP' + A Pa 



aPP' - AP A 



5 TP r + AP A 



^aPP' In 



aRR' - AP A 



aPP' + AP A 



(237) 



We now integrate expression <|237|) in r, in the interval from r = to r = oo. The first term on the right hand side 
of Eq. ([237!) vanishes due to the falloff conditions (see Eqs. lfT94ll - i[202|) and Eqs. (f204|) - <[212]l ^ . We then obtain the 
following expression 



/•OO /"OO 

/ dr {Pa&A + P B SR + P r ST) dr (P M SM + P R SR + P Q SQ) = Suj [A, P, V, P A , P r ] , 

Jo Jo 



(238) 



where 6u> [A, P, T, Pa, Pr] is a well defined functional, which is also an exact form. This equality shows that the 
difference between the Liouville form of {P, A, T; Pr, Pa, Pr} and the Liouville form of {R, M, Q; P r , Pm, Pq} is 
an exact form, which implies that the transformation of variables given by the set of equations l|236j) is canonical. 

Armed with the certainty of the canonicity of the new variables, we can write the asymptotic form of the canonical 
variables and of the metric function F(t, r). These are, for r — > 

F{t,r) = 4R 2 (t)A {ty 2 r 2 + 0(r 4 ) , (239) 

R(t,r) = R (t) + R 2 (t)r 2 + O(r 4 ), (240) 

M{t,r) = M a {t) + M 2 (t)r 2 + 0(r 4 ) , (241) 

Q(t,r) = Q {t)+Q 2 (t)r 2 + 0(r 4 ) , (242) 

PR(t,r) = 0{r), (243) 

P M (t,r) = 0{r), (244) 

P Q (t,r) = 0(r), (245) 



1 



with 

M (t) 

M 2 (t) 
For r — > oo, we have 



a- 1 P - 1 (<)g2(t) + -a 3 Pg(t) 
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(246) 



^aR (t)R 2 {t) (3a 2 R Q {t) 



64R 2 (t)A (ty 2 ) +a- x Q {t)Ro \t) Q 2 (t) - -Q (t)R^(t)R 2 (t) .(247) 



F(t,r) = — r 2 -2(7 1 (t) + 2p(t))r- 1 +0°°(r- 2 ). 

R{t,r) = r + 16p(t)a- 2 r- 2 + O co {r- 3 ), 

M{t,r) = M + (t) + 0°°(r- 1 ), 

Q{t,r) = Q+itj + O^ir- 1 ), 

P R (t,r) = 0°°(r- 4 ), 

P M (t,r) = 0°°(r- 6 ), 

P Q (t,r) = 0°°(r- 2 ), 



(248) 

(249) 
(250) 
(251) 
(252) 
(253) 
(254) 
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where M + (t) = a(r](t) + 3p(i)), as seen before in the surface terms (see Eq. (|214p ). 

We are now almost ready to write the action with the new canonical variables. It is now necessary to determine 
the new Lagrange multipliers. In order to write the new constraints with the new Lagrange multipliers, we can use 
the identity given by the space derivative of M, 

M' = -A- 1 (R'H + a- l R- l P K (H r - TG) + a - l AR- l P T G) . (255) 

Solving for H and making use of the inverse transformations of A and Pa, in Eqs. (|234|) and l|235p . we get 

H = _M'F-^ + FP M P R + a -^F-^QG 
(F-i(R') 2 -F^ 5 

H r = P M M' + P R R' + P Q Q' , (257) 
G = -Q' . (258) 

Following Kuchaf [Tfl|, the new set of constraints, totally equivalent to the old set H(t,r) = 0, H r (t,r) = 0, and 
G(t,r) — outside the horizon, is M'{t,r) — 0, P R (t,r) — 0, and Q'(t,r) — 0. By continuity, this also applies on the 
horizon, where F(t, r) = 0. So we can say that the equivalence is valid everywhere. So, the new Hamiltonian, the 
total sum of the constraints, can now be written as 

NH + N r H r + $G = N M M' + N R P R + N Q Q' . (259) 

In order to determine the new Lagrange multipliers, one has to write the left hand side of the previous equation, Eq. 
(|259p . and replace the constraints on that side by their expressions as functions of the new canonical coordinates, 
spelt out in Eqs. (|256p - (|258[) . After manipulation, one gets 

N M = r + N rp Mj (260) 

(F-^Rr-FPiy 

nR = NFPm _ + N r Rlj (261) 



(F-i(R') 2 -FPf f , 



1 

NQ = a -iNRTiF-iXQ i+ ^ 

(F-^rty-FPiy 

Using the inverse transformations Eqs. I|234p - (|235p . and the identity R = R, we can write the new multipliers as 
functions of the old variables 

N M = -NF^R'A- 1 + a- 1 N r F- 1 R- 1 AP A , (263) 
iV R = —a~ 1 NR~ 1 P\ + N r R' , (264) 
N Q = a- 1 NA- 1 R- 1 F- 1 R'Q~N r T-a' 2 N r R- 2 F- 1 AP A P r -^>, (265) 

allowing us to determine its asymptotic conditions from the original conditions given above. These transformations 
are non-singular for r > 0. As before, for r — > 0, 

N M {t,r) = -iWiWAoWifeCt)- 1 + 0(r 2 ) , (266) 
N R (t,r) = 0(r 2 ), (267) 
N Q (t,r) = ia- 1 iV 1 (t)Qo(i)Ao(t)i? ( 7 1 Wi?^ 1 W-$o + 0(r 2 ), (268) 



and for r^oowe have 



N M (t,r) = -A>+(t)+0°°(r- 4 ), (269) 
N R {t,r) = 0°°(r- 2 ), (270) 
N Q (t,r) = -$+(t) + oo (r- 1 ). (271) 
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These conditions (|266p - l|27ip show that the transformations in Eqs. (|263p - (|264j) are satisfactory in the case of r — > oo, 
but not for r — » 0. This is due to fact that in order to fix the Lagrange multipliers for r — > oo, as we are free to do, 
we fix N+(t), which we already do when adding the surface term 



dtN+M 4 



(272) 



to the action, in order to obtain the equations of motion in the bulk, without surface terms. The same is true for 
However, at r = 0, we see that fixing the multiplier N M to values independent of the canonical variables is not the 
same as fixing A^Aq 1 to values independent of the canonical variables. The same is true of the fixation of with 
respect to $o- We need to rewrite the multipliers N M and for the asymptotic regime r — > without affecting 
their behavior for r — > oo. In order to proceed we have to make one assumption, which is that the expression given in 
asymptotic condition of M(t, r), as r — ► 0, for the term of order zero, M = ^a~ 1 Qn(t) 2 R ~ 1 (t) + -^a 3 Ro(t) 3 , defines 
Ro as a function of M and Qo, where Ro is the horizon radius function, Rn = i?j 1 (Af , Qo). Also, we assume that 

3 

Mq > Af crit (Qo) = ^t-- With these assumptions, we are working in the domain of the classical solutions. We can 

34 

immediately obtain that the variation of Rq is given in relation to the variations of M and of Qo as 

-l 



6M - 



QqSQq 
aRo 



(273) 



This expression will be used when we derive the equations of motion from the new action. We now define the new 
multipliers N M and N Q as 



N M = -N 



NQ 



where g(r) = 1 + 0(r 2 ) for r 



M 



3(T 3 

II-.'/) </«."„) (^-iS 



Ql 

2aR 2 Q 



N M gQ 
and g(r) 



3a 3 
"32" 

0° 



multipliers N M and AfQ, have as their properties for r 

N M {t,r) = 



2aR 2 

~ 5 ) for r 

— > 00 



(274) 



(275) 



00. These new multipliers, functions of the old 



N + (t) 



o°° 

O 00 



(r~ 5 ) 



(276) 
(277) 



and for 







N M (t,r) 
N M (t,r) 



N a M (t)+0(r 2 ) 
$o(*) + 0(r 2 ) , 



(278) 
(279) 



where Nq is given by 



N, 



A I 



1 



1 



When the constraints M ' = 



2 u z V 2 

Q' hold, the last expression is 



32 



(280) 



ATM = NlA ~l 



(281) 



With this new constraint A^ M , fixing A^Aq 1 at r = or fixing N^ 1 is equivalent. The same happens with fixing $0 

or at r — * 0. There are no problems with N^, which is left as determined in Eq. I|26ip . The new action is now 
written as the sum of Ss, the bulk action, and Sgs, the surface action, 



M, R, T, P M ,P R , Pr;N M , N R , N Q 



= / dt 



M 



dr {p M Af + P R R + P Q Q + N Q Q' - N R P R 



N 



(l-g)M'+g 



3a^ 
16 



Ql 

2aRl 



(aR M' - QoQ') 



+ 



dt ( ^aR 2 N™ - N+M+ + $ Qo - 



(282) 
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The new equations of motion are now 



M 


= o, 


(283) 


R 


= N^, 


(284) 


Q 


= o, 


(285) 


Pm 


- (^v M )', 


(286) 


P R 


= o, 


(287) 


Pq 


= (N Q y, 


(288) 


M' 


= o, 


(289) 


P R 


= o, 


(290) 


Q' 


= 0. 


(291) 



Here we understood N M to be a function of the new constraint N M , defined through Eq. I|274p , and to be a 
function of the new constraint , defined through Eq. (|275p . The resulting boundary terms of the variation of this 
new action, Eq. <|282[1 . are, first, terms proportional to 5M, SR, and 5Q on the initial and final hypersurfaces, and, 
second, the term 

J dt (J^aRlSN* 1 - M+SN+ + Q 5$ - Q+SQ+^J (292) 

Here we have used the expression in Eq. (|273|1 . The action in Eq. I|282p yields the equations of motion, Eqs. (|283[i - 
(|29ip , provided that we fix the initial and final values of the new canonical variables and that we also fix the values 
of Nq , N+, $o, and $+. Thanks to the redefinition of both the Lagrange multipliers, from N M to N M and from 
to , the fixation of those quantities, Nq' 1 , N+, $o, and <f>+ has the same meaning it had before the canonical 
transformations and the redefinition of N M and . This same meaning is guaranteed through the use of our gauge 
freedom to choose the multipliers, and at the same time not fixing the boundary variations independently of the choice 
of Lagrange multipliers, which in turn allow us to have a well defined variational principle for the action. 



D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion 
(|283p ~ l|29ip allow us to write M and Q as independent functions of the radial coordinate r, 



M{t,r) = m(t), 
Q(t,r) = q(t). 

The reduced action, with the constraints and Eqs. (|2Q3|i - <|294[) taken into account, is 



m, p m . q, p q ; iV M , N+, $ , $ 



= / dt 



p m m + p q q h, 



where 



Pm 



Pq 



dr P M , 



drPc 



and the reduced Hamiltonian, h, is now written as 

1 



h(m, q;i) = --aR^Nf + N+m + I $+ - $ J q 



(293) 
(294) 



(295) 

(296) 
(297) 

(298) 



with i? n being the horizon radius. We also have that m > Af crit (q) = according to the assumptions made in 
the previous subsection. Thanks to the functions Nq (t), N + (t), &o(t), and $+(t) the Hamiltonian h is an explicitly 
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time dependent function. The variational principle associated with the reduced action, Eq. (|295|1 . will fix the value of 
m and q on the initial and final hypersurfaces, or in the spirit of the classical analytical mechanics, the Hamiltonian 
principle fixes the initial and final values of the canonical coordinate. The equations of motion are 



m = 0, (299) 

q = 0, (300) 

Pm = «^h " a^j) ^ ~ ^+ ' (301) 

Pq = - q (|« 3 4-^) + (302) 



The equation of motion for m, Eq. <|299[1 . is understood as saying that m is, on a classical solution, equal to the mass 
parameter M of the solutions in Eq. I|177p . The same holds for q in Eq. I|300p . i.e., q is the charge parameter in Eqs. 
(|177p ~ l|179p . In order to interpret equation of motion for p m , Eq. I|30ip . we have to recall that from Eq. (|231|) one 
has Pm = —T', where T is the Killing time. This, together with the definition of p m , given in Eq. I|296p . yields 

Pm = To T+ , (303) 

where To is the value of the Killing time at the left end of the hypersurface of a certain t, and T + is the Killing time 
at spatial infinity, the right end of the same hypersurface of t. As the hypersurface evolves in the spacetime of the 
black hole solutions, the right hand side of Eq. <|300[1 is equal to T — T + . Finally, after the definition 

Pq = Co - £+ , (304) 

obtained from to Eqs. (|217j> . <|233[1 . and (|297[1 . Eq. ([302P gives in the right hand side £ — £+, which is the diference 
of the time derivatives of the electromagnetic gauge £(t, r) at r = and at infinity. 



E. Quantum theory and partition function 



Following Sec. IIIII we need not repeat here the steps taken there in order to build the time evolution operator. 
Thus the operator K can now be written as 



K (m, q; T, 8, 3 , 5+) = cxp 



-imT + ^ai? n 6 - iq (S+ - S ) 



(305) 



The same composition law as in Sec. IIIII holds here as well. Here again T is the Killing time elapsed at right spatial 
infinity and 9 is the boost parameter elapsed at the bifurcation circle; So and S + are line integrals along timelike 
curves of constant r and constant angular variables, at r = and at infinity. 



F. Thermodynamics 

We can now build the partition function for this system. The path to follow is to continue the operator to imaginary 
time and take the trace over a complete orthogonal basis. Our classical thermodynamic situation consists of a three- 
dimensional spherically symmetric charged black hole, asymptotically anti-de Sitter, in thermal equilibrium with a 
bath of Hawking radiation. Ignoring back reaction from the radiation, the geometry is described by the solutions in 
Eq. (|177p and (j 1 7Q[) . Thus, we consider a thermodynamic ensemble in which the temperature, or more appropriately 
here, the inverse temperature (3 is fixed. This characterizes a grand canonical ensemble, and the partition function 
Z{f3, (j>) arises naturally in such an ensemble. To analytically continue the Lorentzian solution we put T = —i(3, and 
6 — 2ni, this latter choice based on the regularity of the classical Euclidean solution. We also choose So = and 
S+ = i{34>. 

We arrive then at the following expression for the partition function 

Z(fl, $) = Tr \K(-if3, -2m, 0, i04>)] . (306) 
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From Eq. Q305P this is realized as 



Z(f3,(f>)= I ^idmdqexp -/3(m - q(f>) + iraR? (m|m) . (307) 



A 



Since (m|m) is equal to 5(0), one has to regularize (|307|1 . Again, following the Louko- Whiting procedure we 
have to regularize and normalize the operator K beforehand. This leads to 



-Zxen(/3, 4>) — N J \i dmdq exp -/3(m - q<£) + TnxRjH , (308) 



where M is a normalization factor and A is the domain of integration. Provided the weight factor fi is slowly varying 
compared to the exponential in Eq. (|308p , and using the fact that the horizon radius is function of m and q, the 
integral in Eq. l|308j) is convergent. Changing integration variables, from m to i?^, where 



1 Jn3 , ^ 



the integral Eq. l|308j) becomes 



m -T2 aR ^^ (309) 



Zren(A <t>) = M I JidR^dq exp(-J.) , (310) 
J A' 



where A! is the new domain of integration after changing variables, and the function ^(i?^, q), a kind of an effective 
action (see is written as 



M/? h .q) := ±0 a 3ij3+AiL_ ; , , q ;rn /,'- . {311) 



The domain of integration, A', is defined by the inequalities < R^ and q 2 < 3a 4 i?^/16. The new weight factor 
ju includes the Jacobian of the change of variables. Since the weight factor is slowly varying, we can estimate the 
integral of Z Ten ((3, <f>) by the saddle point approximation. The critical points are given as a pair of values (i?^, q) 



K = w^V 1 ^^ 2 ^' (312) 

q* = aR±$. (313) 

The critical points belong to the domain of the effective action (|31 ip , and so the effective action has a null derivative at 
these points in Eqs. (|31 2j) - <|3 13[) . As 4> 2 > -^/3 _2 a _2 0~ 2 7r 2 there are always two critical points, given by (|313p - (|312|) . 
The critical point at (RT , q~) is outside the domain of physical interest as R7 < 0. The critical point at (R^: q + ) is 
an global minimum. The expression l|313p implies that the critical points are on the straight line defined by itself, on 
the plane of the arguments of the effective action (|31ip . Then, if we replace the value of the charge as a function of 
i. e., q(-Rjj) as given in l|313p into the effective action we obtain a one variable function of R^, written as 

7 * ( R h) = ^ a3R h ~ l^ R h ~ a * R i > ( 314 ) 

which is a function of R^ on the straight line defined by q = aR^4>. The zeros of this effective action on the straight 
line are given by 

4 = 0, (315) 



R^ = \&p~ 1 a- 2 n\l±\jl + ^l3 2 a 2 ft^j . (316) 

Here R^ < is outside the physical domain of interest. As (ft 2 > — yg/3~ 2 a~ 2 0~ 2 7r 2 there are always three zeros of 
the function <|314[) . which means that there are no phase transitions, as the zeros of the effective action maintain their 
respective hierarchy whatever the value of <f>, which is the order parameter of the system. So the classical solution 
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is found in the global minimum of the critical points (-R^,q + ) of the effective action i*(i? n ,q) (Eq. I|31ip ). The 
partition function Z Ten (/3, <j>) can now be written, through the saddle-point approximation 

^ren(/3,0) « Pexp [-/*(#+,q+)] , (317) 
lnZ ren (/3,0) « In P q+) . (318) 

where P is a slowly varying prefactor, and j q + ) m the effective action evaluated at the global minimum critical 

point. By ignoring the prefactor's logarithm, which closer to (-R^,q + ) is less relevant, we are able to determine the 
value of m at the critical point, where we find that it corresponds to the value of the mass of the classical solution of 
the black hole given in Eq. l|177p . Thus, when the critical point dominates the partition function, we have that the 
mean energy (E) is given by 

(E) = ~ lnZ ren « ^(B+) a " \ a ^K = m+ ' ( 319 ) 
where m + is obtained from Eq. I|309p evaluated at (R-^, q + )- The thermal expectation value of the charge is 

(Q) =/r 1 i=lnZ re n«q + , (320) 

0(f> 

as given in Eq. (|313|1 . We now can write the temperature T = /3 _1 as a function of the critical point given by R^ 
and q + 

T = ^-(saRt-^ictyUR+Y 3 ) , (321) 



4?r V h 16 

with i?^ explicitly written in l|312p . and where q + is given in Eq. (|313p . If we choose the maximum charge of 

the domain given by q 2 < 3a 4 i?^/16, then we find ourselves with the extreme solution of the black hole, whose 

temperature, if q 2 = 3a 4 i?^/16 is replaced into (|321[) . is equal to zero. It can be shown that dm + /d/3 < 0, which 

through the constant <j> heat capacity = —(3 2 {d{E) /d/3) tells us that the system is thermodynamically stable. 
Finally there is the entropy of the black hole, S, which is given by 

S = (l - P-^j (In Z ren ) « 4™i?+ 2 . (322) 

Here we see that despite the fact that the extreme solution has null temperature, it has an entropy given in (|322p . 
which is independent of the value of the charge in the domain of validity of the approximation. 

V. HAMILTONIAN THERMODYNAMICS OF A REPRESENTATIVE CHARGED DILATONIC BLACK 

HOLE (U = -3) 

A. The metric 

For lo = — 3, the general metric in Eq. ([2]), the (f) field in Eq. ©, and the vector potential iJH) reduce to the 
following, respectively, 



ds 2 



(aRf - 2(aR)i M - ^{aR)Q 2 



dT 2 + ^ +R 2 d<p 2 , (323) 

(aR) 2 -2{aR)iM- \{aR)Q 2 



e" 2 ^ = —r, (324) 
A = --(aR)iQdT, (325) 
with 2M — b and a — 3^|A| = ^^-l^ 1 . In Fig. [3] the Carter-Penrose diagram of the black hole solution for the case 
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R=0 




R=0 

FIG. 3: The Carter-Penrose Diagram for the lu — — 3 case, with M > 0. 



lu = — 3 is shown, which is analogous to the case u) = 0, where again the singularity at R = is a curvature singularity, 
and R = oo is spatial infinity. The horizon radius is given by the largest real solution of the following equation 

(a R h ) 2 - 2 (a R h )$ M - i(a i? h )Q 2 = 0. (326) 

Note that this diagram, Fig. [31 holds for static black holes with positive mass M > 0. Here I and I' are the right 
and left static outer regions of the Carter-Penrose diagram, respectively, and II and IF are the future and past inner 
dynamical regions of the Carter-Penrose diagram. i? n is the horizon radius given in Eq. I|326p and the spacelike 
infinity is asymptotically anti-de Sitter, both in the right and in the left regions, denoted by I and F, respectively. 
There are no extreme black hole solutions for lu = — 3 [251 ]. 

Since we are now familiar with the whole formalism, we will be briefer in this section omitting several of the details. 



B. Canonical formalism 



The action becomes, following Eqs. lf2l]) - lf2"3|) . and up to surface terms, 



S[A, R, A, R; N, N r ] = 



dt 



dr \ 4(aRy^NX 2 AR- (a R)~i N^AR + -{a R)~ 2 N^AR^R 2 



+ {aR)-^N~ 1 R(N r A)' + (a R)~ ? N' 1 N r AR 1 - (a R)~* N' 1 N r (N r )' AR' 
-{aR)-^N' 1 (N r ) 2 A'R' - ^(a R)-? N' 1 N r AR' 1 RR' + ^(a R)~? N-\N r ) 2 AR- 1 (R') 2 

+ ^{aR)-^NA- 1 R- 1 {R') 2 - {a R)~* N^A' 1 )' R' - (aRy^NA^R" 



2N- 1 A- 1 {aRy^R[t - 



(327) 



Depending on the situation we use three different letters containing the same information but with slightly different 
numerical values. Thus, a, A, and I are related by a 
the action above, we obtain the conjugate momenta 



numerical values. Thus, a, A, and I are related by a 2 — -yA 2 = |i 2 , where I is defined now as the AdS length. From 



P A = -N- L (aR)-2 [R- R'N 



1 



P R = -N- L (aR)-2 i A- (AN r y - -AR- L (R- R'N r )} , 



P r = 4A^ 1 A- 1 (ai?)-2 j R(r-$' 



(328) 
(329) 
(330) 



By performing a Legendre transformation we obtain the Hamiltonian, which is a sum of constraints, i.e., 



H = A \-{aR)ip A P R + ^aiR-iAPl + (aRY 



(A-^'R' + A^R" - ^A-^R'fR- 1 
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A{aR)-?\ 2 AR+ -NAR^iaR^Pf 
+ N r {P R R' - AP' A - TP r } + $ {-P r } = NH + N r H r + $G . (331) 
We can now write the action in Hamiltonian form, which reads 

S[A, R, T, P A , Pr, Pr; N, N r , $] = J dt dr jp A A + P R R + P r t - NH - N r H r - <Ig} , (332) 

with the constraints defined in Eq. (|331[) . Prom here we derive the equations of motion for the canonical variables 
and respective canonical momenta 



A 
R 
f 

Pa 



-N(aR)*P R + -NaiR-2AP A + (N r A)' . 
-N{aR)?P A + N r R' , 
-NiaR^AR^Pr + (N r P r )' + , 



NaiR-^Pl 



(N(aR)-^ R' 



A 2 ~ ^N(a R)-^ (R'j 2 R^A-' 2 + AN(a R)-i X 2 R + N r P A 



-N(a R)^ R~ 1 P 2 , 



(333) 
(334) 
(335) 

(336) 



)-Na?R-?P K P R + iNa^R-iAPl - Na~i (A -1 )' (R~i)' 
2 8 



]N(aR)-3A- 1 R- 1 R" - ( (NaT*)' A^BT^) - ^Nar^A- 1 ^^)" 



Pr = 



+2Na-?\ 2 AR-? + (N r P R )' + —NA(aR)^R- 2 P 2 , 
N r P T . 



(337) 
(338) 



For a correct variational principle to be applied, we have to find out what surface terms are left over from the variation 
performed, with the purpose of deriving the equations of motion in Eqs. P33p ~ p38p . These surface terms are, 



Surface terms = N(a R)~iR'A~ 2 6A - N(a R)~i A^SR' + |(aii)^j' A^SR + ^(a R)~? A^R^R'SR 



- N r P R SR + N r A6P A + N r T6P r + $6P r 
In order to know the form of (|339[i for r — ► 0, we assume 



A(i,r 

R(t,r 
PA(t,r 
P R (t,r 

N(t,r 
N r (t,r 

T(t,r 
P r (t,r 



A + O(r 2 ), 

Ro + R 2 r 2 + 0(r 4 ) , 

0(r 3 ) 

0(r), 

iVi(t)r + 0(r 3 ) 

0(r 3 ), 

0(r), 

Qo + Q 2 r 2 + 0(r 4 ) , 
*oW + 0(r 2 ). 



The surface terms become 



Surface terms| r=0 = -Ni(aR) ?A 1 SR - <&o<>Qo 



The asymptotic conditions for r — > oo are assumed as 

V6 



A(t,r) 



Ir- 1 + 2?l2r)(t)r-2 + 0°°(r- 3 ) 



(339) 



(340) 
(341) 
(342) 
(343) 
(344) 
(345) 
(346) 
(347) 
(348) 



(349) 



(350) 
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R(t,r) 




Pa (t r) 


- 0°°(Y~ 2> ) 






N(t, r) 




N r (t, r) 


— 0°°(r~ 2 ) 

— \i ) , 


T(t, r) 


= 0°°(r- 2 ), 


Pr(t,r) 


= Q+W + O^r- 1 ), 


$(t,r) 





where I is the AdS length. The surface terms, Eq. I|339p . for r — + oo, are written as 

Surface termsl^^ = N+5M+ + $+<5Q+ , 

where 

M+{t) = 2 5 3~3 rj(t) + 2 2 33 p(t) . 



(351) 
(352) 
(353) 
(354) 
(355) 
(356) 
(357) 
(358) 



(359) 



(360) 



Here we pay attention to the fact that we can define the above quantity M + , in Eq. (|360p , because we do as explained 
in subsection III Al when uj < — 1 and absorb the divergence in the definition of M. This means that the M + above 
in Eq. I|360p is in fact related to the mass plus the divergence, defined above as M(i2*), see Eq. ff20|) and below in 
subsection IV CI The surface term added to (|332|) is then 



Sdi 



A,R,Q ,Q + ;N,$o,$ + = / dt [2 a - 1 (a Rq)? 1 - N+M+ + $ Qo - ®+Q+ 



(361) 



With this surface term added we obtain a well defined variational principle. What remains after variation of the total 
action, Eq. ([332]) and Eq. pel} , is 



dt 



(2a- 1 (ai?o)^(A r iA,7 1 ) - M+5N+ + Q S$ - Q+S$ H 



(362) 



The surface terms coming from the variation of the total action disappear as the result of the fixation of NiA 1 on 
the horizon, r = 0, and of N+(t) at infinity, r — > oo, with the same happening respectively to $o at r = and to $ + 
at r — > oo, 



C. Reconstruction, canonical transformation, and action 



Repeating the steps of the two previous sections, we now give the main results concerning ui = —3. The vector 
potential is 



A = -~(aR)?QdT + £, 



where there is the function £ = £(t, r). Expanding for (t, r) we find 

A = 

out of which we infer that 



^(aR)?QT' + A dr+ (-^-( a R)^QT + i ) ,11 . 



Pr = 



The metric function is given by the expression 



F= (aR) 2 - 2{aR)2 M — -(aR)P 2 
The Killing time T = T(t, r) is a function of (t, r), where we find that 

— T' = P" 1 AP A (aP)^ , 



(363) 



(364) 



(365) 



(366) 



(367) 
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which is equal to minus the conjugate momentum of the new variable M, or Pm 
help of Eq. l|367p we find F = F(t, r) as a function of the canonical variables, 



P M (t, r) 



F = {B!fK- 2 - aRPl 



-T'(t,r). With the 
(368) 



Summing up, the canonical transformations are 

R = R, 
1 

2 ( 
Pr, 



M = -(aR)-i [ (aR) z - -{aR)P 2 - F 
Q 



Pm 

Pq 



RR'AP' A + -(R') 2 P A A - RR"AP A 



Pr-F^R^A 

+ -R- 1 AP A - 4A 2 P~ 1 P A PA + -clF^APaPt , 
4 8 

F~ 1 AP A (aR)i , 
-r+i(aJ?) AP A P r . 



RR'A'P A 



(369) 



In addition, the relevant inverse transformations back to the old variable A and respective conjugate momentum Pa 
are 



A = (F-\R'f-FPlY , 
P A = {aR)-^FP M (F- l {R'f -FP 2 M y h . 
We have to show that this transformation is canonical. This requires using the identity 



(370) 
(371) 



P A SA + P R 8R + P r 8T - P M 8M - P R 8R - P Q 8Q 



-(aR)-*6R In 



(aR)-2R' + AP A 



(aR)-?R'- AP A 



1 



S TP r + AP A + -(aR)-^R'hx 



(aR)-?R' - AP A 



(aR)-2R' + AP A 



|(372) 



We now integrate Eq. I|372p . in r, from r = to r — oo. The first term on the right hand side of Eq. (|372p vanishes 
due to the falloff conditions, Eqs. (|340|) - (|348j) and Eqs. I|350p - (|358p . We obtain then the following expression 



dr (P A SA 



P R SR + P r ST) dr (P M SM + P R SR + P Q 5Q) = 8lo [A, P, T, P A , P r ] , 

Jo 



(373) 



where 5ui [A, P, T, P A , Pr] is a well defined functional, which is also an exact form. As above, this equality shows that 
the difference between the Liouville form of {P, A, T; Pr, P a , Pr} and the Liouville form of {R, M, Q; P r , Pm, Pq} 
is an exact form, which implies that the set of transformations Q369P is canonical. 

With this result, we write the asymptotic conditions of the new canonical coordinates for r — * 

(374) 
(375) 
(376) 
(377) 
(378) 
(379) 
(380) 

where we have 



F(t,r) 


= ARlA^r 2 + 


0(r 4 ), 


R(t,r) 


= R a + R 2 r 2 - 


(-0(r 4 ) 


M(t,r) 


= M Q + M 2 r 2 


+ 0(r 4 


Q(t,r) 


= Qo+Q2T 2j , 


-0(r 4 ) 


P R(t>r) 


= 0(r), 




P M (t,r) 


= 0(r), 






= 0(r). 





M = -(aP ) 



(aP ) 2 --(aPo)Qo 



M 



J2 



\{aRv)Ql + (aP ) 2 ] P 2 2a 2 R 2 n R% - \ [2(aR )QoQ 2 + aQ 2 R 2 ] - 4P 2 A 



4P (aP )^ 



2(aP )^ 



(381) 
(382) 
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and a = I 1 a/8/3. For r — > oo , we have 

4,_„ 2 16 / ^r 1 



F(t,r) = li-V-^p^- (8^+3^)) rU0»(r»), (383) 



R(t, r) 


-oh 


■2 

Pit) 


Q(*,r) 


= g + (<) + o°°(r- 




M(f,r) 


= M + (t)+0°°(r- 






= 0°°(r- 4 ), 






= 0°°(r-i), 






= 0°°{r- 2 ). 





}OQ /„— 1 \ 



(384) 

(385) 
(386) 
(387) 

(388) 
(389) 



where M+(t) is defined in Eq. (|360| . 

We now write the future constraint A/' as a function of the older constraints 



M' = -A _1 [r'H +{aR)^P A (H r -TG)-^{aR)^KP T G S j . (390) 

Using the inverse transformations of A and Pa in Eqs. I|370p and (|371[) . we obtain the same form for the old constraints 
as functions of the new variables 

H = M'F-W + FPmPr - ^iaR^QQ'F^R! 

H r = P M M' + P R R' + PqQ' , (392) 
G = -Q' . (393) 

The new Hamiltonian, the total sum of the constraints, can now be written as 

NH + N r H r + $G = N M M' + N R P R + N Q Q' . (394) 

The new multipliers are, using Eqs. (|391[1 - (|394[) . 

7\/ Tp — 1 T>l 

N M = -+N r P M , (395) 

{F-W-FPlY 

N R = ^ r + , (396) 

(F-i(R') a -^) s 

ArG 1 N{aR)?F- l R'Q Arrn - /onw , 

A/ g = -^-r + N r P Q - $ . (397) 

2 (F-HR') 2 -FPf f ) 5 

Using the inverse transformations Eqs. I|370p - (|37ip . and the identity R = R, we can write the new multipliers as 
functions of the old variables 

N M = -NF^R'k- 1 + N r F^ 1 KP^a R)^ , (398) 
N R = —NP\(a R)^ + N r R' , (399) 

N Q = — —N(aR)^F~ 1 R'A^ 1 Pr — N r T + ^N r (aR) AP^Pr — $ . (400) 



For r->flwe have, 



N M (t,r) = --N^AoRz 1 + 0(r 2 ), (401) 



iV R (i,r) = 0(r 4 ), (402) 
W Q (i,r) = -iQ ^ 1 (aPo)^iA -$o + 0(r 2 ). (403) 
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and for r-»oowe have 



N M (t,r) = -N + {t) + 0°°{r-'), 
N R (t 7 r) = 0°°(r-i), 
N Q {t,r) = -#+(<) + O co (r- 1 ). 



(404) 
(405) 
(406) 



Again, for r — > 0, fixing N M (t, r), which means fixing NiA R 2 , is not equivalent to fixing A^Aq . It is thus necessary 
to rewrite N M for r — > 0. The same is true of fixing N®(t, r) when r — > 0. In this case it is not the same if one fixes 
— jQqR^ 1 {a Rq)i Nik — $o or just $o- So, assuming Mo as a function of i?o and Qo allows one to define the horizon 
radius Rq = R^(A1q,Qo). We are thus working in the domain where Mq > M ct[t (Qo) = (a 3 — 1) 12i yg , the domain 
of the classical black hole solution. The variation of Rq is given in terms of the variations of Mq and of Qo in the 
expression 



'3 i 1 

SR = ( -a(aR)? - -a(aR y?Q 2 



The new multiplier N is obtained from the old N as 



SM + -(aRo^QoSQo ) ■ 



(407) 



N M = -N 



M 



[ I - ,/ i + f a /?„ ) -• <i [ ja(aRo)>-^a (a R ) I Q 2 



-i -i 



where a — ^p|A| = I 1 ^/8/3, and g(r) = 1 + 0(r 2 ) for r — * and g(r) — 0°°(r 5 ) for r — ► oo. This new multiplier, 
function of the old multiplier, N M , has as its properties for r — > oo 



1 ,~ 
2 



(408) 
(409) 



and as its properties for r — > 



AT M (i,r) 
A> Q (i,r) 



iV M (t,r) 
iV M (t,r) 



A> + (t) + 0°°(r- 2 ), 
l> + (t) + oo (r- 1 ), 



iV M (i) + O(r 2 ), 
$ (i) + O(r 2 ), 



where A/g is given by 



N, 



A I 



A^iA i? 2 

M'=Q=Q' 



3 3 

— a 2 Rn a Qn 

8 48 ^° 



(410) 
(411) 



(412) 
(413) 



(414) 
(415) 



With this new constraint N M , fixing AiA 1 or fixing Nq 1 is equivalent, and the same happens with N®, where fixing 
$o is fixing the multiplier in the limit r — > 0. There are no problems with 7V^, which is left as determined in Eq. 



(13991 . 

The new action is then, summing both the bulk and the surface terms, 



S 



M, R, Q, P M ,P R , Pq; N m ,N R , N q 



dt 



-N 



A I 



dr {P M M + P R R + PqQ + N Q Q' - A^ R P R 



(oi2o)-»Af' + QoQ / l]} + 

(2a- 1 (ai? )2A> M - A+M+ + $ Qo - 



(416) 
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The new equations of motion are now 



M 


— n 

— u , 




it 


atR 


(418) 


Q 


= o, 


(419) 


Pm 


- {N M y, 


(420) 


P R 


= o, 


(421) 


PQ 


= (N^y, 


(422) 


M' 


= o, 


(423) 


P R 


= o, 


(424) 


Q' 


= 0. 


(425) 



Here we understood N M to be a function of the new constraint, defined through Eq. I|408p . The resulting boundary 
terms of the variation of this new action, Eq. l|416p . are, first, terms proportional to SM, SR and 5Q on the initial and 

final hypersurfaces, and, second, the term J dt ^2a _1 (a Ro)i SNq' 1 — M + dN + + QoS&o — Q+<5$+^ . Here we used the 

expression in Eq. (@07j). The action in Eq. (|416p yields the equations of motion, Eqs. (|41T|l - (|425fl . provided that 
we fix the initial and final values of the new canonical variables and that we also fix both the values of and 
N+, and $o and Thanks to the redefinition of the Lagrange multiplier, from N M to N M , and from to 
the fixation of those quantities, Nq 4 and N + , and $o and has the same meaning it had before the canonical 
transformations and the redefinition of N M and . This same meaning is guaranteed through the use of our gauge 
freedom to choose the multipliers, and at the same time not fixing the boundary variations independently of the choice 
of Lagrange multipliers, which in turn allow us to have a well defined variational principle for the action. 



D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion, 
Eqs. (|417p - (|425p . allow us to write M and Q as independent functions of the radial coordinate r 



M(t,r) = m(t), 
Q(t,r) = q(t). 

The reduced action, with the constraints and Eq (|427p taken into account, is 



m, q, p m , p q ; N* 1 , N + ,$ , $ 



= / dt 



p m m + p q q h, 



where 



Pm 



Pq 



drP M 



drP c 



and the reduced Hamiltonian, h, is now written as 

h(m, q; t) = -2a' 1 (a R h ) i 7V A/ + iV+m + (l>+ - * ) q , 

with i?^ being the horizon radius. We also have that m > M crit (q). The equations of motion are then 

m = 0, 
q = 0, 

'3 



Pm 



(aR h Y 



■a(aR h )* - -a(aR h ) 



p q = ^i\2 a ( aR h) 2 ~ z a ( aR h) 2 q^ 



(426) 
(427) 



(428) 

(429) 
(430) 

(431) 

(432) 
(433) 

(434) 
(435) 
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Here m and q are equal to the mass parameter M and the charge parameter Q of the classical solution in Eq. (|323p . 
The third equation of motion, Eq. I|434p . describes the time evolution of the difference of the Killing times on the 
horizon and at infinity, due to p m = T — T + and Eq. I|429p . The last equation, Eq. <|435|i . gives the time evolution of 
Pq as the difference between the time derivatives of the electromagnetic gauges £ (i, r) on the horizon and at infinity 
due to p m = £o — £+ and Eq. (|430p . 



E. Quantum theory and partition function 



The steps developed in Sections HUl and [TV] can be readily used here. So we do not spell out the corresponding 
calculations in detail. Thus, obtained as above, the time evolution operator is 



K (m, q; T, 0, So, S+) = exp —imT + 2ia 1 (a i? n ) 2 — iq (S+ — 
We can now develop the thermodynamic analysis as in Sections IIIII and IIVI 



(436) 



F. Thermodynamics 

We can now build the partition function for this system, with the Hamiltonian given in Eq. (|431[) . The steps are 
the same as was the case with u> = oo and u = 0. The thermodynamic ensemble is also the grand canonical ensemble. 
Again, with T = —if), = — 2ni, S = 0, and S + = i(3(f>, we write the general form of the partition function as 



Z = Tr 

This is realized as 



K(-i/3,-2m,O,0(p) . (437) 



ZveniP, 4>) = N f jLdR h dq exp(-i*) , (438) 

J A' 

where N is a normalization factor, A! is the domain of integration given by < and q 2 < 60^, and the function 
I* C^h' a kind of an effective action [lfj]. q), is written as 



/.(i2 h ,q) := ^(oJZ h )* - r|o*q + 47ro-*)i^-q0/3. (439) 



Its critical points are at 

i?± = ^tt/T^- 2 |l± ^l + ^(3 2 a 2 n- 2 4> 2 j , (440) 
q± = -2^(ai?±)"i. (441) 

Given that (f> 2 > — |/3 _2 a _2 7r 2 there are always two critical points, but only Rt is physical as it is positive, Rt > 0. 
The other critical point is negative, R^ < 0. Determining the zeros of the action <|439[1 when Eq. <|441[) holds, we get 



R^ 1 = Anp- 1 ^ 2 M± ^1 - i/3 2 a 2 7r- 2 2 j . (442) 

Now, for different values of 4> 2 we may have different signs for the action <|439[1 evaluated at the physical critical point 
i?^. So, for (f> 2 > 8f3~ 2 a~ 2 n 2 there are no zeros and I*(Rt) > 0, where we take q + = q + (i?^) to be a function of the 

critical point radius Rt. If 4> 2 — 8/3~ 2 a~ 2 7r 2 then the effective action <|439|) has one zero only, at Rt, i. e., I*(R^) = 0. 

Finally, if (f> 2 < 8/3~ 2 a~ 2 7r 2 then we have two zeros, given in Eq. (J442J) , and I*(i?^) < 0. Then, despite the variation 

of the sign of the effective action at the physical critical point (ii^,q + ), there are no new minima when the order 

parameter <j> is changed across the border given by <j) 2 — 8f3~ 2 a~ 2 Tr 2 , which means that there are no phase transitions, 
as the classical phase corresponds always to the global minimum. We can now derive the basic thermodynamic results, 
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as long as the approximation for the saddle point is valid, which means that we have to work in the neighborhood of 
the classical solution (|323p . where the critical point dominates. The expected value of the energy E is 



IE) [ -— + 13-^ 1 lnZ ren « 2 (ai V % " ^ aR V^ = m ' ^ 



The charge expectation value is 

(0) = /3- 1 T n lnZ ren «q + . (444) 

Here we have the temperature T = given by 

1 f$ 



= i{r 2R i-i«^r)- («») 

It can be shown that the derivative of m + (/3) with respect to /? is negative, which means that the heat capacity is 
positive, i. e., — — (3 2 ( ^ > 0. The system is thus stable. Finally the entropy is given by 



d 



S=[l- P T ) (In Zreni/3)) « 4tt y/a^R+ . (446) 
The result recovers the entropy for the three-dimensional charged dilatonic black hole with ui — — 3 (see also [28l|). 



VI. CONCLUSIONS 



We have extended the work done in [l[ by adding a Maxwell term to the three-dimensional Brans-Dicke action, and 
thus, have continued Louko and collaborators' program @, OJ 0, H, S 0] (see also [H,@|)) of studying the thermodynamic 
properties through Hamiltonian methods of black holes in several theories in different dimensions. Specifically we 
have calculated the thermodynamic properties of black hole solutions with asymptotics that allow a well formulated 
Hamiltonian formalism in three-dimensional Brans-Dicke-Maxwell gravity. Only certain values of the Brans-Dicke 
parameter w are allowed in this juncture. The corresponding theories are general relativity, i.e., ui — > ±oo, a dimen- 
sionally reduced cylindrical charged four-dimensional general relativity theory, i.e., uj = 0, and a theory representing 
a class of charged dilaton-gravity theories, with a typical w given by w = —3, all of these coupled to a Maxwell field. 
We have shown that for these theories, with well defined asymptotics for geometric and electromagnetic variables, 
the formalism fits well. It was possible to have a well defined electric field for all the values of uj studied. As in [l| 
several modifications were needed. First, in the powers of the fall-off conditions, and second due to the presence of 
the scalar field, which was reflected in the fact that it changed the powers of the radial R coordinate, to name a few. 
Of course, the addition of electromagnetism also introduced a new pair of canonical variables, which added a good 
deal of complexity to the calculations. Namely, the treatment of the mass term when electromagnetic interaction is 
included is not trivial at all in this Hamiltonian formalism in three dimensions, since there is a divergent asymptotic 
behavior of the mass M. Within the Hamiltonian formalism we solved this problem by introducing a renormalized 
mass. This redefinition of the mass fits perfectly into the formalism for the u — » ±oo (i.e., general relativity), it is 
not needed in the uj — case, and conforms to the uj = — 3 case. Although, in order to build a three-dimensional 
Lorentzian Hamiltonian theory with electric charge, these modifications had to be made, a quantum theory and a 
statistical description of the systems in question in a grand canonical ensemble were derived, and the corresponding 
thermodynamics, with precise values for the entropy of the black holes studied, was established. 
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